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I 

I 

I  SECTION  I 

INTRODtCTION 

Currently,  the  North  Warning  System  (NWS)  employs  the  Analog  Devices  Corporation 
(ADC)  12-bit,  4-MHz  analog-to-digital  (A/D)  convener  model  MOD- 1205.  During  the  full- 
scale  development  of  the  AN/FPS-124  Radar,  the  ADC  convener  was  selected  as  the  only 
one  that  met  the  system  performance  requirements  (six  of  them  arc  used  in  the  three  receiver 
channels  of  each  radar).  This  convener  has  been  in  operation  for  three  years  with  no  known 
problems;  however,  Paramax  (formerly  Unisys  Corporation)  is  proposing  to  replace  it  w  ith 
model  ADC-001 15-600,  manufactured  by  Data  Device  Corporation  (DDC).  Paramax 
delivered  two  ADC  conveners  (referred  to  as  Units  1  and  2)  and  a  DDC  convener  for  our 
evaluation  (shown  in  figures  1(a)  and  1(b)),  enabling  us  to  comment  on  their  engineering 
change  proposal. 

As  pan  of  our  assessment  and  evaluation  of  the  DDC  A/D  convener,  we  conducted  a  scries 
of  tests  to  verify  the  Unattended  Radar  (UAR)  receiver  B2  performance  requirements  1 1 1  for 
the  baseline  ADC  unit  (idcndfied  in  table  1 ).  A  performance  comparison  revealed  that 
ADCs  Unit  1  was  extremely  noisy  and  did  not  meet  spurious-free  dynamic  range  (SFDR) 
radar  B2  specifications.  The  dynamic  test  on  Unit  2  revealed  that  it  was  noisier  than  Unit  1 
and  did  not  meet  radar  B2  Specifications  cither  (note  that  both  have  been  in  operation  for 
three  years).  Then,  using  the  ADC  tests  as  a  baseline,  the  relative  performance  of  the  DDC 
unit  was  assessed  and  we  also  verified  if  it  met  the  B2  requirements.  This  report  presents  t 
detailed  analysis  of  the  evaluation  and  dynamic  performance  testing  of  these  three 
conveners. 


1.1  A/D  CONVERTER  TESTS 


The  dynamic  performance  (nonlinearities  and  noise)  of  an  A/D  convener  can  be 
characterized  by  a  spectral  analysis  test  (2,3,4).  consisting  of  a  series  of  sinusoidal  signals  of 
different  frequencies  and  amplitudes.  The  input  frequency  selected  for 


Table  1.  UAR  Receiver  B2  Performance  Requirements 


Voltage  Range  ±  2.048  Volts 

Conversion  Time  Data  valid  275  +  25  nanoseconds  after  third  encode 

command 


Integral  Linearity  ±0.1  percent  of  full-scale  ±  0.5  least  significant  bit  (LSB)* 

Monotonicity  At  1.25  MHz  -  no  missing  codes,  each  code  edge-to-edgc 

width  >tcn  percent,  <190  percent  of  nominal 


Alternating  Current  At  0.54  MHz  -  spurious  70.0  dB  below  full-scale 

(ac)  Linearity 


Signal-to-Noise  Ratio  At  0.54  MHz,  63.0  dB  minimum  root  mean  square  (RMS) 

(SNR)  signal-to-RMS  noise 


Transient  Response  Transient  response  (for  full-scale  step  input):  attains  full- 

scale  accuracy  to  ±  5.0  mVolts  in  200  nanoseconds  maximum 


Overvoltage  Recovery  Time  Overvoltage  recovery  time  (for  2x  full-scale  input):  recovers 

+  5.0  mV  (millivolts)  of  12-bit  accuracy  within 
200  nanoseconds 


♦  NOTE:  One  LSB  is  defined  as  the  resolution  of  an  A/D  converter,  which  is  equal  to 

(2*^  X  the  full-scale  voltage  of  the  converter),  where  N  =  the  number  of  bits  of 
the  converter. 


the  test  is  representative  of  the  expected  signal/clutter  frequencies.  For  our  test,  we  select 
seven  frequencies,  two  of  which  are  specified  in  the  radar  B2  specification  [1].  Full- 
amplitude,  single-tone,  sine-wave  testing  is  the  most  stressful  test  for  A/D  converters  [3] 
because  the  full-amplitude  signal  level  maximizes  the  slew  rate  of  the  input,  and  tests  all 
encoding  levels.  The  results  of  the  spectral  analysis  test  are  used  to  determine  the  effective 
nun^r  of  bits,  the  SFDR,  and  the  harmonic  distortion  of  the  A/D  converter. 

The  B2  specification  requires  that  the  A/D  converter  be  tested  for  monotonicity  at  a  test 
frequency  of  1.25  MHz,  and  that  the  separation  of  all  of  the  code  edges  between  adjacent 
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codes  be  1.0  ±  0.9  least  significant  bits  (LSB),  which  is  a  specification  of  differential 
linearity.  It  is  difficult  to  measure  monotonicity  or  code-edge  separation  with  a  high- 
firequency  test  signal;  either  special  test  equipment  or  elaborate  data  processing  is  required. 
TTie  monotonicity  test  is  usually  done  at  direct  current  (dc).  which  docs  not  reflect  high- 
firequency  performance.  Although  it  is  a  variant  of  the  sine-wave  test,  we  employ  custom 
software  to  perform  monotonicity  analysis  of  data  obtained  with  a  pure  sine-wave  input 
voltage. 

The  transient  and  overvoltage  recovery  performances  of  the  converters  are  determined  by 
applying  a  pulsed  waveform  as  input  and  determining  the  time  required  to  recover  from  the 
pulse.  Ideal  pulses  have  no  overshoot  or  ringing,  and  negligible  rise  and  fall  times.  Since 
commercial  test  equipment  used  in  our  test  facility  is  unable  to  produce  a  pulse  which  is 
clean  to  12  bits  accuracy  with  a  400  ohm  load,  a  special  test  fixture  was  designed  and  built 
(the  details  are  described  in  section  5). 

A  systematic  approach  to  characterizing  high-speed  A/D  converters  based  on  sinusoidal 
testing  was  developed  as  part  of  a  MITRE-sponsored  research  project  [4].  To  test  the  A/D 
converters  digitally,  a  laboratory  test  facility  was  constructed.  The  next  section  presents  the 
key  parameters  used  to  characterize  them  dynamically. 

Section  2  of  this  repon  briefly  describes  the  DDC  and  ADC  converters  that  are  tested,  the 
test  equipment,  and  the  data  analysis  algorithm  used  to  characterize  them.  Section  3  presents 
the  dynamic  performance  results  of  the  DDC  converter  (tested  at  seven  frequencies)  and  the 
two  ADC  converters.  Section  4  provides  the  details  of  the  monotonicity  performance,  while 
section  5  presents  the  transient  and  overvoltage  recovery  test  results.  Appendix  A  presents 
the  general  sine-wave  estimation  algorithm;  appendix  B  is  a  FORTRAN  listing  of  the 
simplified  sine-wave  estimation  algorithm;  and  appendix  C  presents  the  statistical  analysis 
which  is  used  in  evaluating  the  monotonicity  test  results. 
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TEST  DEFINITION 


This  section  briefly  defines  the  dynamic  performance  measures  [5,6]  used  to  characterize  the 
A/D  converters.  The  quantitative  values  of  the  parameters  are  determined  in  the  frequency 
domain  by  performing  spectral  analysis. 

12.\  Signal4o>Noise  Ratio 

A  standard  measure  of  A/D  converter  performance  for  signal  processing  applications  is  the 
ratio  of  the  signal  power  to  the  noise  power  in  the  Nyquist  frequency  band  used,  usually 
referred  to  as  the  signal-to-noise  ratio  (SNR)*.  The  SNR  is  defined  as 

Sm  =  10  \og{PslNs )  in  dB,  ( 1 ) 

where  /j  is  the  signal  power,  and  is  the  total  noise  power  in  the  Nyquist  frequency  band, 
which  extends  from  the  low  frequency  to  half  the  sampling  frequency. 


1.2.2  Signal-to<Noise*and*Distortion  Ratio 


Signal-to-noise-and-distortion  ratio  (SINAD)  is  defined  as  the  ratio  of  the  signal  power  to  the 
total  error  power.  The  total  error  power  includes  the  power  in  the  harmonics  as  well  as  the 
power  caused  by  random  noise. 


SINAD  =  10 


log 


in  dB, 


(2) 


where  Pg  is  the  signal  power,  N^  is  the  noise  power,  and  P^  is  the  power  in  the  harmonics. 


1  Nyquist  bandwidth  is  . 
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1.2.3  Effective  Number  of  Bits 


The  theoretical  dynamic  range  for  an  ideal  A/D  convener  with  a  full-scale  sine-wave  input  is 

6.02  N  +  1.76  dB  [5,6],  where  N  is  the  number  of  bits  (quantization  noise  is  the  only  noise). 
For  a  nonideal  A/D  converter,  the  SINAD  and  the  effective  number  of  bits,  are  directly 

related  by  the  following  equation; 


5/NAD  (in  dB)-1.76 
6.02 


(3) 


The  effective  number  of  bits  is  computed  from  the  SINAD  measured  in  the  frequency 

domain  using  spectral  analysis.  The  number  of  bits  lost  is  the  difference  between  the 
nominal  bits  of  the  A/D  converter  and  the  effective  number  of  bits.  Note  that  JVgjf  is  a 

function  of  the  input  test  signal  frequency  and  amplitude;  it  is  merely  a  different  way  of 
expressing  the  SINAD. 

1J1.4  Spurious-Free  Dynamic  Range 

Spurious-free  dynamic  range  (SFDR)  [4]  is  defined  as  the  ratio  of  the  power  of  the  sine- 
wave  input  to  the  power  of  the  harmonic  with  the  highest  amplitude,  measured  in  dB,  up  to 
half  the  sampling  frequency  and  excluding  the  dc  component. 

U.5  Integral  Nonlinearity 

The  integral  nonlinearity,  as  shown  in  figure  2,  is  the  maximum  deviation  of  the  A/D  transfer 
function  from  a  best-fit  straight  line. 

12.6  Differential  Nonlinearity 

The  differential  nonlinearity,  also  shown  in  figure  2,  is  a  measure  of  how  much  the  code 
widths  vary  from  the  ideal  code  width  of  1.0  LSB. 
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Binary  Output  Code  Binary  Output  Code 


1L2672-1 


Figure  2.  Enors  in  A/D  Converters 
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\2.1  Code  Edge 


We  define  the  code  edge  [7]  of  output  code  i  as  the  lowest  input  voltage  V ,  so  that  there  is  at 
least  a  50  percent  probability  that  the  output  code  equals  or  exceeds  i. 

12.S  Monotonicity 

An  A/D  converter  is  monotonic  [7]  if  the  expected  mean  output  code  never  decreases  for 
increasing  input  voltage,  and  if  it  never  increases  for  decreasing  input  voltage.  It  is  strictly 
monotonic  if,  for  any  finite  increase  (or  decrease)  of  the  input,  the  expected  mean  output  also 
increases  (or  decreases).  Generally,  a  monotonic  A/D  convener  is  also  strictly  monotonic 
unless  it  has  missing  codes  which  cause  a  finite  region  of  input  voltage  to  produce  the  same 
output.  The  B2  specification  requires  that  the  converters  for  the  NWS  be  strictly  monotonic. 


12  SUMMARY  AND  CONCLUSIONS 

As  pointed  out  in  the  introduction,  we  obtained  a  new  DDC  A/D  convener  and  two  three* 
year-old  operational  ADC  A/D  conveners  from  Paramax  which  were  subjected  to  an 
extensive  series  of  dynamic  tests.  The  objective  of  this  study  was  to  assess  the  DDC  A/D 
converter  for  the  NWS.  Using  the  ADC  tests  as  a  baseline,  the  relative  performance  of  the 
DDC  unit  was  assessed  and  we  also  verified  if  it  met  the  UAR  receiver  B2  performance 
requirements. 

This  report  briefly  describes  the  characteristics  of  the  two  A/D  converters  tested  and  the 
digital  test  bed  used  to  characterize  them.  We  analyze  the  data  in  the  frequency  and  phase 
domains.  Based  on  the  spectral  analysis  test,  we  evaluate  these  key  A/D  converter  dynamic 
parameters;  the  noise,  effective  number  of  bits,  SFDR,  SNR,  and  the  SNDR  for  seven 
different  frequencies  and  eight  different  amplitude  levels  at  each  frequency.  The 
monotonicity  of  the  converters  is  tested  at  a  frequency  of  1.25  MHz  [1].  The  monotonicity 
test  is  sensitive  to  converter  noise,  so  to  reduce  the  chances  of  drawing  erroneous 
conclusions,  we  employ  three  different  methods  of  analysis  (histogram  method,  code-mean 
method,  and  code-edge  method)  to  determine  the  differential  nonlinearity  and  monotonicity 
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performance  of  the  converters.  The  transient  and  overvoltage  recovery  tests  are  performed 
with  a  custom-designed  pulse  generator  circuit.  As  a  result  of  our  investigation,  we  reach  the 
following  conclusions: 

•  The  DDC  converter  exhibits  its  best  noise  performance  (0.7  to  0.8  LSB  RMS 
noise)  at  the  B2  specified  frequency  band  from  0.54  to  1.25  MHz.  For  full- 
amplitude,  very  low-frequency  (i.e.,  100  kHz),  and  high-frequency  (i.e., 

1.5  MHz)  signals,  the  noise  of  the  converter  is  0.95  LSB  and  1.0  LSB, 
respectively.  For  signals  near  the  Nyquist  frequency  of  2.0  MHz,  the  noise  is 
higher  at  all  signal  amplitude  levels. 

•  The  DDC  converter  loses  approximately  1.3  to  1.5  bits  with  respect  to  an  ideal 
12-bit  converter,  which  corresponds  to  10.5  to  10.7  effective  bits  under  most 
dynamic  conditions.  However,  for  full-scale  amplitude  signals,  it  loses 
approximately  2.2  bits  at  both  the  high  (>1.5  MHz)  and  low  frequency 

(0.1  MHz). 

•  The  SFDR  of  this  converter  is  71.0  to  73.0  dB  for  large  amplitude  signals.  At 
other  frequencies  it  is  similar  to  the  results  at  0,54  and  1.25  MHz,  except  at  very 
low  frequencies  of  1(X)  kHz  and  270  kHz,  where  the  distortion  is  higher  and  the 
SFDR  is  65.0  dB  for  large  signal  amplitudes. 

•  The  integral  nonlinearity  of  the  DDC  conveners  is  (under  ±  1 .0  LSB) 
exceptionally  good  at  a  nominal  frequency  of  544  kHz.  At  two  lower  frequencies 
of  101  and  271  kHz  the  integral  nonlinearity  is  less  than  ±  1.5  LSB,  which  is 
three  times  better  than  the  B2  specification. 

•  Since  many  subranging  A/D  converters  behave  differently  with  positive-going  or 
negative-going  signals,  we  present  the  monotonicity  analysis  for  positive-going, 
negative-going,  and  combined  slopes  to  gain  better  insight.  The  DDC  converter 
is  strictly  monotonic  and  marginally  meets  the  radar  B2  specification 
requirements;  it  appears  to  have  several  codes  which  may  be  slightly  wider. 
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Confirming  this  with  high  confidence  would  require  more  samples,  greater  than 
the  512  K*  used  in  our  analysis. 

•  The  DDC  converter  does  not  meet  the  transient  and  overvoltage  recovery  B2 
performance  requirements. 

•  The  two  operational  three-year-old  ADC  converter  units  do  not  meet  most  of  the 
radar  B2  specifications;  however,  during  qualification  testing  of  the  radar  three 
years  ago,  they  did  meet  them. 

•  Unit  1  has  1.15  to  1.4  LSB  of  RMS  noise  for  full-scale  input  signals,  which  is  an 
excess  of  0.6  LSB  of  noise  in  comparison  to  the  DDC  converter.  This  converter 
loses  2.4  to  2.9  bits  for  full-scale  input  signals,  and  2.6  to  3.0  bits  for  signals  that 
are  30  dB  below  full-scale.  The  SFDR  of  this  converter  is  60  to  65  dB  for  large 
signals. 

•  Unit  2  has  1.1  LSB  of  RMS  noise  at  most  of  the  frequencies  (0.27  MHz, 

0.54  MHz,  and  1.5  MHz)  and  for  different  signal  amplitudes  at  each  frequency. 
However,  at  a  signal  frequency  of  0.1  MHz,  the  noise  rises  appreciably  to 
3.0  LSB  which  may  be  associated  with  the  missing  code  problem.  TTtis  converter 
loses  between  2.2  and  2.5  bits  at  all  frequencies  and  different  signal  amplitudes  at 
each  level.  However,  for  a  full-scale  amplitude  signal  with  a  frequency  of  1 .5 
MHz,  it  loses  3.5  bits,  which  may  be  again  due  to  the  missing  code  problem.  The 
SFDR  of  this  converter  is  62  to  66  dB  for  large  signals. 

•  The  spectrum  plots  of  both  units  show  an  increase  of  6  dB  in  noise  floor  over  a 
wideband  in  the  upper  half  of  the  Nyquist  band,  possibly  caused  by  high-order 
nonlinearity  which  results  from  second-stage  error  correction  circuitry  not 
performing  properly. 


*  throughout  this  report  1  K  =  1024  samples. 
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•  The  ADC  converter  is  nonmonotonic,  it  has  gross  errors  (particularly  for 
negative-going  voltages),  and  it  does  not  meet  the  radar  B2  specification. 

•  The  ADC  converter  does  not  meet  the  overvoltage  and  transient  recovery 
performance  B2  requirements  which  were  not  a  part  of  the  B  2  specification 
during  the  qualification  testing  of  the  radar. 

•  The  results  of  the  Paramax/DDC  and  MITRE  evaluations  of  the  DDC  A/D 
converters  are  summarized  in  table  2.  The  DDC  converter  meets  all  B2 
performance  requirements  and  exceeds  the  performance  of  the  two  three-year-old 
ADC  units,  with  the  exception  of  the  overvoltage  and  transient  recovery  time, 
which  is  out-of-specification.  This  is  not  expected  to  present  a  problem,  since  the 
A/D  converter  in  the  system  (receiver)  is  preceded  by  a  hard  limiter  that  is  set 
slightly  above  the  maximum  full-scale  input  level.  This  will  prevent  the  A/D 
converter  from  being  driven  to  twice  die  full-scale  input.  Since  the  transient 
recovery  time  is  only  slightly  out-of-specification,  and  the  settling  characteristic 
of  this  converter  has  a  smooth  exponential  decay,  the  only  system  impact  is  a 
moderate  attenuation  of  very  high-frequency  signal/clutter  components. 

13.1  System  Impact 

In  this  section  we  present  a  preliminary  analysis  [8,9,10]  to  evaluate  the  effect  of  converter 
noise  on  the  system.  The  available  dynamic  range  depends  on  the  sensitivity  of  the  system  to 
harmonics,  the  noise  produced  by  the  A/D  converter  with  typical  input  signals,  and  the 
system  noise  prior  to  conversion.  For  illustrative  purposes,  we  calculate  a  noise  figure,  that 
is,  the  degradation  of  the  SNR  expressed  in  dB  caused  by  imperfections  of  the  A/D  converter. 
Table  3  shows  the  effect  of  the  two  A/D  converters  with  receiver  noise  set  at  1.5  quanta.  As 
shown,  the  DDC  converter  normally  produces  an  additional  0.7  dB  due  to  the  internal  noise 
over  and  above  the  loss  from  ideal  quantization  noise  (0. 16  dB).  However,  the  three-year-old 
operational  ADC  A/D  converter  produces  an  additional  loss  of  1.0  dB  compared  to  the  DDC 
converter. 
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Table  2.  DDC  Model  ADC  0015-600  A/D  Converter  Performance  Summary 
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Since  the  radar  is  in  operation  with  the  ADC  converter,  and  appears  to  be  fully  functional 
despite  its  converters  being  out-of-specification,  we  recommend  an  additional  study  to 
evaluate  the  system  impact  of  not  meeting  the  specifications. 

Hie  new  DDC  converter  meets  most  the  B2  Specifications,  as  did  the  ADC  converter  units 
during  qualification  testing  of  the  radar.  Our  test  results  revealed  that  the  three-year-old 
ADC  converter  units  did  not  meet  any  of  the  specifications.  We  recommend  that  the  DDC 
converter  performance  be  checked  over  an  extended  period  of  time,  after  being  in  operation 
for  roughly  three  years,  to  ensure  that  it  will  not  impact  the  system  performance. 

The  present  test  methods  for  the  convener  based  on  sine  waves  are  not  adequate  because  it  is 
difficult  to  correlate  the  results  with  real-world  radar  performance.  Perhaps  new  test  methods 
employing  simulated  radar  signals  may  be  required  to  quantify  the  degradation  effects 
clearly,  and  provide  additional  insight  showing  how  a  system's  parameters  will  dictate 
specifications  for  a  converter. 


Table  3.  System  SNR  Loss  Caused  by  Two  A/D  Converters 


Converter  Type 


WQ 

Analog  Devices 

Receiver  Noise  (in  LSB) 

1.5 

1.5 

A/D  Converter  Noise  (in  LSB) 

0.7-1.3* 

1. 1-3.0** 

Total  Noise  (in  LSB) 

1.65-1.98 

1.86-3.35 

Loss  (in  dB) 

0,86-2.43 

1.87-6.98 

♦  Noise  is  less  than  0.8  LSB  RMS,  except  with  high-frequency  (>1.5  MHz)  signals 
within  S.O  dB  of  saturation.  This  is  not  a  normal  input  condition. 

*♦  The  3.0  LSB  figure  is  produced  by  the  errors  caused  by  missing  codes,  which  result 
from  high-frequency  (>1.0  MHz)  signals  with  an  amplitude  near  full-scale,  which  is 
not  a  normal  input  condition. 
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SECTION  2 


DESCRIPTION  OF  TV/0  A/D  CONVERTERS  AND  THE  A/D 
INSTRUMENTATION  FACILITY 

2.1  CHARACTERISTICS  OF  TWO  NWS  A/D  CONVERTERS 

Cuirendy,  the  NWS  uses  the  ADC  12-bit,  5-MHz  A/D  converter  model  MOD- 1205, 
employing  a  two-stage  digital  correcting  subranging  architecture  with  pipelining.  The  digital 
output  has  a  delay  of  three  clock  cycles,  plus  0.275  microseconds  (psec)  after  the  input  is 
sampled.  At  a  4-MHz  sampling  rate,  the  digital  output  has  a  total  delay  of  1.025  psec.  The 
converter  has  a  nominal  linear  range  from  -2.048  to  +  2.047  volts  (V),  and  a  resolution  of 
12  bits  (which  corresponds  to  0.(X)1  volt  (one  millivolt)).  The  input  impedance  is  selectable 
by  an  external  connection  to  be  either  50.0  or  400.0  ohms. 

DDC  manufactured  model  ADC-001 15-600  as  a  replacement  (illustrated  in  figure  1(b)). 
Because  it  is  intended  to  be  a  replacement,  it  is  based  on  a  printed  circuit  board  which  has  the 
same  size,  mounting  provisions,  connector  pinout,  and  location  as  the  ADC  unit.  However,  it 
is  based  on  a  small  hybrid  circuit  assembly  which  is  a  complete  12-bit  A/D  converter.  In 
addition  to  the  hybrid  assembly,  a  few  components  produce  an  input/output  interface  that  is 
identical  to  the  ADC  converter.  Much  of  the  DDC  converter  circuit  board  is  empty  in 
comparison  to  the  ADC  convener.  Table  4  shows  the  primary  specifications  of  both. 


2.2  DESCRIPTION  OF  THE  DIGITAL  TEST  BED 

To  study  A/D  converter  limitations  on  signal  processing  algorithms  used  in  communication 
and  radar  systems,  we  developed  a  systematic  approach  to  characterizing  high-bandwidth 
and/or  high-resolution  A/D  conveners  [4].  The  investigation  analyzed  the  technical  issues 
involved  in  accurately  evaluating  them  in  the  frequer  cy  and  phase  domains  [7],  including  the 
construction  of  a  versatile,  high-performance,  A/D  instrumentation  facility. 
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Table  4.  Summary  of  the  Most  Important  Characteristics  of  the  A/I)  Converters* 


Manufacturer 

ADC 

D12C 

Bits 

12 

12 

Conversion  Rate  (MHz) 

5.0 

5.0 

Full-Scale  Range  (±2.098  Volts) 

4.096 

4.096 

Data  Output  Pipeline  Delay  (Cycles) 

3.0 

3.0 

Data  Output  Delay  (psec) 

0.275 

0,275 

Input  Impedance  (Ohms) 

50/400 

50/4(X) 

Overvoltage  Range  (Volts) 

±  4.0 

±4.0 

Supply  Voltages 

±  15, +  5, -6 

±  15, 5.  -6 

Total  Power  (W) 

13.0 

10.0 

SNR  (dB)  for  0.54  MHz  Signal 

66.0 

64.0 

SFDR  (dB) 

from  dc  to  1.0  MHz 

70.0 

70.0 

1.0  MHz  to  2.5  MHz 

65.0 

65.0 

Transient  Recovery  Time  (psec) 

0.2 

0.2 

Settling  Accuracy  (Millivolt) 

0.5 

0.5 

*Ba$ed  on  manufacturer's  specification 


Digital  testing  requires  taking  blocks  of  digital  data  from  the  A/D  converter,  storing  them  in 
real  time  in  a  data  ^quisition  memory,  and  then  transferring  them  to  a  computer,  where 
various  analyses  are  performed  in  temporal,  spectral,  and  statistical  domains.  This  prevents 
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any  reconstruction  errors  that  can  occur  when  the  digital  samples  are  converted  back  to 
analog  values  via  a  digital-to-analog  (D/A)  converter. 

Figure  3  illustrates  digital  testing  of  an  A/D  converter  in  block  diagram  form.  A  detailed 
description  of  the  A/D  instrumentation  capability  at  MITRE-Bedford  is  listed  in  table  5.  The 
core  of  the  instrumentation  facility  consists  of  very  high-speed,  high-capacity  (dual-input  or 
single-interleaved  mode)  data  acquisition  memories.  One  data  acquisition  memory  can  log 
data  records  up  to  two  megawords  (Mwords)  in  a  dual  channel  mode  or  four  Mwords  in  an 
interleaved  mode.  A  second  data  acquisition  memory  is  expandable  up  to  16  (duji)  or  32 
(interleaved  mode)  Mwords.  Each  input  channel  can  accept  data  words  up  to  16  bits  wide  at 
a  rate  up  to  100  Mwords  per  second  (200  in  interleaved  mode).  The  two-channel  capability 
is  a  convenient  feature  for  the  characterization  of  A/D  interfaces  such  as  digital  radios  and 
radar  systems  with  complex  inphase  and  quadrature  outputs.  The  data  acquisition  system  is 
supported  by  two  personal  computers  (PCs)  for  data  analysis  and  software  development.  The 
analog  sources  shown  in  figure  3  are  just  a  sample  of  the  type  of  signal  generators  available; 
they  represent  the  kind  of  analog  input  signals  that  the  A/D  converter  might  encounter  in  a 
real-world  communication  or  radar  environment  (e.g.,  narrowband  and/or  broadband,  signal, 
jammer,  noise,  or  clutter).  Their  outputs  are  further  conditioned  for  better  noise  and/or 
linearity  with  passive  filters  before  being  applied  to  the  A/D  unit  under  test.  A  photograph  of 
the  A/D  instrumentation  laboratory  is  shown  in  figure  4. 

The  sine-wave  is  applied  to  the  A/D  converter  under  test  and  the  resulting  digital  output  is 
stored  in  the  Tektronix  High-Speed  Data  Cache,  a  large,  fast,  sequentially-addressed 
memory.  A  Compaq  PC  controls  it  via  an  IEEE-488  bus,  employing  DADiSP,  a  general- 
purpose  signal  processing  program.  A  macro  command  entered  at  the  keyboard  commands 
the  Data  Cache  to  store  64  K  (65536)  consecutive  data  points,  and  then  transfer  them  to  the 
Compaq  PC.  The  data  record  is  displayed  on  the  screen,  and  then  stored  as  a  file  on  a  floppy 
disk. 
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Figure  3.  A/D  Instrumentation  Facility  Generic  Block  Diagram 
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Table  5.  A/D  Instrumentation  Capability 


Function 
Signal  Generator 


A/D  Instrumentation  Hardware  and  Software  Description 

HP-3335A(2).  HP-3325A.  TF-2096.  HP-8662A.  HP-8644A. 
HP8656B,  KH4400A. ... 


Signal  Conditioner  Passive  bandpass  and  low-pass  filters  (30,  HF-1  high-intercept 

two-tone  combiner) 


Data  Acquisition 


Data  Processing 


Software 


Tek  9503, 100-Msps  (ECL),  25-Msps  (TTL),  16-bii  (two-channel), 
2-Mword/channel 

Tek  9504, 100-Msps  (ECL),  25-Msps  (TTL),  16-bit  (two-channel), 
1 6-Mword/channei 

MITRE  custom  memory,  50-Msps  (ECL),  16-bit  (one-channel), 
4-kword/channel 

Compact  Deskpro  386/33  MHz,  200-Mb  hard  disk,  16-Mb  RAM 

Compact  Deskpro  386S/25  MHz,  80-Mb  hard  Disk,  16-Mb  RAM 
HP-9386  (2),  Laser  Printer 

Digital  Signal  Processing  Spreadsheet  (DADiSP) 

Excel,  Turbo  Pascal,  Turbo  C++,  Axum  Graphics,  IEEE488, ... 
MITRE  custom  A/D  characterization  software  and  macros 


2.3  PHASE-PLANE  ANALYSIS  OF  CONVERTERS  WITH  SINE  WAVES 

In  many  signal  processing  applications,  the  dynamic  performance  of  an  A/D  converter  is 
characterized  by  the  spectral  analysis  test,  which  consists  of  applying  a  series  of  sinusoidal 
signals  of  different  frequencies  and  different  amplitudes,  and  analyzing  the  frequency 
spectrum  of  the  output  to  characterize  the  nonlinearities  and  noise.  It  does  give  accurate 
measurements  for  the  large-scale  statistics,  such  as  harmonic  distortion  and  total  noise  of  the 
A/D  converter,  but  it  is  difficult  to  detemine  the  error  sources  within  the  device. 
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Phase-plane  analysis  [11]  evaluates  the  errors  and  noise  at  each  code  level  of  the  A/D 
converter  and  provides  insight  into  the  internal  behavior  of  the  device.  It  is  a  way  of 
manipulating  a  mass  of  noisy  data  so  that  measurements  made  under  similar  input  conditions 
can  be  correlated.  By  evaluating  statistical  propenies  of  small  subgroups  of  the  total  data, 
underlying  causes  of  correlations  are  revealed.  This  technique  is  extremely  useful  with  high- 
resolution  A/D  converters.  A  sine  wave  is  applied  to  the  input,  and  the  output  data  is 
recorded.  Using  all  of  these  points,  we  create  an  estimator  that  evaluates  the  exact 
instantaneous  input  voltage  associated  with  each  of  the  measured  points. 

23.1  Sine>Wave  Estimation 

Given  samples  of  a  waveform  which  is  dominated  by  a  sine  wave  with  small  addidonal 
harmonic  distortion  and  noise,  several  algorithms  can  estimate  the  fundamental  sine-wave 
component.  Peetz  et  al.  [12]  present  a  very  brief  description  of  a  general  algorithm  for 
sine-wave  curve  fitting,  requiring  an  initial  estimate  of  the  signal  frequency  (accurate  to 
five  or  six  significant  digits).  Frothing  et  al.  [13]  simplify  the  algorithm  for  the  special 
case  where  the  signal  frequency  is  known  precisely.  This  requires  phase-locked  oscillators 
for  the  test  setup,  which  may  not  be  available.  Appendix  A  presents  a  general  iterative 
algorithm  which  is  an  improved  version  of  Peetz’s  method,  computing  the  sine  wave 
entirely  from  the  measured  data.  In  our  test  setup,  we  employ  extremely  stable,  phase- 
locked  oscillators  for  the  signal  and  sampling  clock.  The  frequencies  of  the  sine  waves  are 
chosen  so  that  there  is  a  prime  integer  number  of  complete  cycles  in  a  data  record, 
allowing  us  to  use  a  fast  algorithm  to  estimate  the  sine  wave  (a  generalization  of  Frothing's 
method). 

We  start  with  a  set  of  N  measurements,  V'-,  for  i  =  1,2,...,  A'  for  a  waveform  consisting  of  a 
large  fundamental  component,  F^,  plus  measurement  errors  and  noise,  e^.  The  fundamental 

is  determined  by  four  parameters:  the  offset,  M;  amplitude.  A;  radian  frequency,  (o;  and 
initial  phase,  6. 


F^  =  M  +  A  cos{o)i  +0)  i  =  1,2,...,  A 


(1) 
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F^  =  M+B  sin(fi}/)  +  C cos(6)/) 


«  =  1.2 . N 


(2) 


where 


e  =  tan-'(%.). 

(3) 

Vi  =  Fi^ei 

II 

(4) 

Let  us  define  K,  K^,  and  as 

N 

(5) 

i=sl 


N 

Ks=  IVi&inicDi)  (6) 

i=l 


N 

Kc  =  '^ViC0s{(0i).  (7) 

i=l 

Substituting  equations  (2)  and  (4)  into  equations  (4)  to  (6),  and  assuming  that  the  e/s  are 
uncorrelated  with  the  sin{ct)i)  and  cos((ui),  results  in  three  simultaneous  equations  in  three 
unknowns:  M,  B,  and  C: 


K  =  Ms  +  Bsg  +  Csc 

(8) 

Ks  —  Msg  +  Bsgg  +  CSjj. 

(9) 

(10) 
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where 


N  N 

j  =  5^sin(Git);  Sc  =  ]^cos(£W) 

i=l  i=l 

N  N  N 

Ssc  =  5^sin(e«)cos(Q>i);  Sgs  =  ^cc  “  ^cos^{Q)i).  (11) 

1=1  1=1  1=1 

The  three  unknowns  M,  B,  and  C  are  determined  by  solving  the  three  simultaneous  equations. 
However,  they  require  that  the  errors,  e,,  must  be  uncorrelated. 

2.3J.1  Spectral  Analysis  Tests 

For  the  spectral  analysis  tests,  the  signal  amplitude  is  chosen  to  maintain  the  A/D  converter 
within  its  linear  range,  and  the  frequency  is  chosen  so  that  there  is  an  exact  integer  number  of 
cycles  within  the  data  record.  Hence,  all  of  the  measured  samples  are  valid,  and  the  s's  can 
be  computed  independently  of  the  measurements. 

s  =  N'j  Sg  =  Sq  =  Sgc  —  0;  Sgg  = 

The  unknowns  Af,  B,  and  C  ate  computed  from 


1  ^ 
1=1 


(13) 


B  =  ^  =  ^'^ViS\n{o>i) 


^ss 


^r=i 


(14) 
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c= 


Vl  cosicoi). 


(15) 


2.3.1.2  Monotonicity  Tests 

For  the  monotonicity  tests  of  the  DDC  convener  (described  in  section  4),  the  signal 
amplitude  is  chosen  to  exercise  all  possible  output  codes,  and  hence  drive  the  A/D  convener 
into  saturation;  the  output  is  a  clipped  sine  wave,  and  the  errors  are  large  at  the  signal  peaks. 
In  our  data  analysis  we  eliminate  the  saturated  samples  from  the  data  record,  so  M,  B,  and  C 
cannot  be  computed  from  equations  (13),  (14),  and  (15).  The  frequency  is  chosen  so  that 
there  is  an  exact  integer  number  of  cycles  within  the  data  record.  The  output  codes  of  the 
unsaturated  samples  have  an  absolute  value  of  <  2047.  Using  all  of  the  unsaturated  samples 
within  the  data  record,  the  coefficients  M,  B,  and  C  are  determined  by  solving  the 
simultaneous  equations  (8),  (9),  and  (10).  The  results  of  the  analysis  are  presented  in 
section  4. 

23.2  niase-Plane  Data  Analysis 

For  each  data  point,  the  input  voltage  and  the  phase  of  the  input  sine  wave  are  calculated. 

The  measurements  and  associated  errors  are  resorted  in  order  of  the  signal  phase  angle  from 
0  to  2n:,  grouping  together  the  measurements  made  under  similar  input  conditions,  but  at 
different  times  [11].  If  multiple  data  records  are  available,  the  data  can  be  combined  into  this 
table. 

Hiis  sorting  maintains  separation  of  the  two  ways  of  having  the  same  input  voltage  with 
positive  and  negative  slopes.  By  averaging  the  error  over  a  small  region  in  phase  space,  we 
obtain  an  estimate  of  the  typical  error  under  the  input  condition  at  the  center  of  that  region. 
This  is  a  measure  of  the  device  distonion  at  that  input  signal  phase.  The  spread  of  the 
measurement  errors  within  a  region,  as  calculated  by  their  standard  deviation,  is  a  measure  of 
the  device  noise.  This  allows  for  the  determination  of  an  estimator  of  the  device  noise  for 
different  signal  phases  and  voltages.  It  is  frequently  possible  to  associate  one  or  more 
regions  of  high  noise  with  specifrc  problems  in  the  device's  internal  construction.  The 
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regions  of  high  noise  with  specific  problems  in  the  device’s  internal  consuuction.  The 
complete  test  requires  a  series  of  measurements  at  different  frequencies  and  amplitudes.  The 
results  presented  in  sections  3  and  4  employ  the  phase-plane  analysis. 
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SECTION  3 


SINGLE-TONE  TEST  RESULTS 


As  discussed  in  section  1.1,  A/D  converter  performance  is  measured  at  the  frequency  of 
interest  by  sampling  a  sine  wave  and  performing  a  fast  Fourier  Transform  (FFT)  on  the 
digital  output  data.  It  is  generally  desirable  to  test  an  A/D  converter  at  many  different  signal 
frequencies  and  at  different  signal  amplitude  levels  at  each  frequency.  As  the  input 
frequency  is  increased,  the  measured  SNR  of  the  A/D  converter  tends  to  decrease.  This 
high-frequency  performance  degradation  of  the  SNR  is  indirectly  related  to  the  signal's  input 
slew  rate. 

Since  the  sampling  rate  of  the  A/D  converter  is  4.0  MHz,  we  test  both  converters  with  an 
input  signal  frequency  near  the  Nyquist  frequency.  We  test  the  DDC  converter  at  seven 
different  frequencies  of  0.1  MHz,  0.27  MHz,  0.54  MHz,  1.0  MHz,  1.25  MHz,  1.5  MHz,  and 
2.0  MHz,  respectively,  and  at  eight  different  amplitude  levels  of  0  dB,  -1  dB,  -3  dB,  -6  dB, 
-12  dB,  -18  dB,  -24  dB,  and  -30  dB,  relative  to  a  full-amplitude  signal  that  toggles  all  the 
bits.  The  sampling  frequency  for  every  case  is  4.0  MHz.  We  test  the  ADC  Unit  1  A/D 
converter  at  two  different  frequencies,  and  we  test  the  ADC  Unit  2  A/D  converter  at  four 
different  frequencies. 

Based  on  the  spectrum  of  the  data,  we  determine  the  dynamic  parameters,  effective  number 
of  bits,  SNR,  and  SFDR,  as  discussed  in  section  1.2.  We  also  compute  the  dynamic 
parameters  of  the  converters  by  generating  a  mathematically  perfect  sine  wave  from  the 
digitized  data,  as  described  in  section  2.  The  difference  between  the  digitized  data  and  the 
perfect  sine  wave  is  the  error  associated  with  each  data  point.  We  then  analyze  the  errors  to 
evaluate  distortion,  noise,  bits  lost,  and  the  dependence  of  the  errors  upon  signal  voltage. 

We  present  plots  of  the  number  of  bits  lost,  SFDR,  noise,  and  SINAD  versus  input  signal 
level,  and  some  typical  FFT  spectral  responses. 
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The  B2  performance  requirements  for  the  baseline  converter  are  identified  in  table  1.  In  this 
section,  we  are  only  validating  the  RMS-signal‘to-RMS  noise  ratio,  SFDR,  and  integral 
linearity  requirements.  The  other  B2  requirements  are  verified  in  sections  4  and  5. 


3.1  TEST  RESULTS  FOR  THE  DDC  12.Brr,  4.MHZ  A/D  CONVERTER 

Our  analysis  of  the  test  results  for  the  DDC  12-bit,  4-MHz  A/D  converter  at  seven  different 
input  signal  frequencies  is  presented  on  the  following  *agcs,  in  the  form  of  plots,  as  follows: 
noise  versus  input  signal  level  under  dynamic  conditioiis  (figure  5);  SINAD  versus  input 
signal  level  (figure  6);  bits  lost  versus  input  signal  level  (figure  7);  SFDR  as  a  function  of 
input  signal  level  (figure  8);  and  the  FFT  spectral  response  for  the  full-amplitude  and  the 
30  dB-below-full-amplitude  sine-wave  signals  (figures  9(a)  through  9(n)). 

Table  6  lists  the  measured  dynamic  performance  at  seven  different  input  signal  frequencies. 
This  converter  exhibits  its  best  noise  performance  at  the  B2  specified  frequency  band  from 
0.54  to  1.25  MHz.  At  these  frequencies,  the  RMS  noise  is  0.7  to  0.8  LSB  for  full-scale  input. 
For  signal  amplitudes  5  dB  below  the  maximum  the  noise  reduces  to  0.6  LSB,  which  is 
7.0  dB  greater  tlian  an  ideal  12-bit  A/D  converter  quantization  noise.  For  a  signal  frequency 
of  0.27  MHz,  the  noise  performance  is  0.7  LSB,  independent  of  signal  amplitude.  For  full- 
amplitude  very-low  frequency  signals  (i.e.,  100  kHz),  the  noise  increases  to  0.95  LSB. 
However,  for  signal  amplitude  10  dB  below  the  maximum  (at  100  KHz)  the  noise 
performance  is  similar  to  that  at  0.54  MHz  (i.e.,  0.6  LSB),  but  at  a  1.5  MHz  signal  frequency, 
it  is  almost  the  same  as  that  of  100  kHz  for  all  amplitude  levels  except  the  full-scale 
amplitude.  For  full-scale  amplitude  signals  at  1.5  MHz,  the  noise  is  slightly  higher  and  is 
equal  to  1.0  LSB.  For  signals  near  the  Nyquist  frequency  of  2.0  MHz,  the  noise  is  higher  at 
all  signal  amplitude  levels.  This  is  approximately  0.75  LSB  for  signals  less  than  15  dB 
below  the  maximum  and  rises  smoothly  to  1.3  LSB  for  signals  at  the  maximum  amplitude. 

Hgure  6  shows  that  the  SINAD  performance  is  extremely  good  at  a  signal  frequency  of 
0.54  MHz,  but  somewhat  poorer  at  1.25  MHz.  For  small  signal  amplitudes  (15  dB  below  the 
maximum)  the  SINAD  is  constant  at  65  dB  for  all  frequencies  and  amplitudes.  However,  for 
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Figure  7.  Bits  Lost  vs.  Input  Signal  Level  for  the  DDC  A/D  Converter 
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Figure  8.  SFDR  vs.  Input  Signal  Level  for  the  DDC  A/D  Converter 
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Figure  9(a).  FFT  Spectral  Response  for  the  DDC  A/D  Converter 
(Input  Signal  Frequency  =  0.10  MHz;  Amplitude  of  the  Signal  =  Full-Scale) 
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Figure  9(b).  FFT  Spectral  Response  for  the  DDC  A/D  Converter 
(Input  Signal  Frequency  =  0.10  MHz;  Amplitude  of  the  Signal  =  -30  dB) 
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Figure  9(c).  FFT  Spectral  Response  fw  the  DDC  A/D  Converter 
(Input  Signal  Frequency  =  0.27  MHz;  Amplitude  of  the  Signal  =  FulFScalc) 
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Figure  9(d).  FFT  Spectral  Response  fw  the  DDC  A/D  (jonverter 
(Input  Signal  Frequency  =  0.27  MHz;  Amplitude  of  the  Signal  =  -30  dB) 
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Figure  9(e).  FFT  Spectral  Response  for  the  DDC  A/D  Convener 
(Input  Signal  Frequency  =  0.54  MHz;  Amplitude  of  the  Signal  =  Full-Scale) 
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Figure  9(f).  FF  l  Spectral  Response  for  the  DDC  A/D  Convener 
(Input  Signal  Frequency  =  0.54  MHz;  Amplitude  of  the  Signal  =  -30  dB) 
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Figure  9(g).  FFT  Spectral  Response  for  the  DDC  A/D  Converter 
(Input  Signal  Frequency  =  1.0  MHz;  Amplitude  of  the  Signal  =  Full-Scale) 
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Figure  9(h).  FFT  Spectral  Response  for  the  DDC  A/D  Converter 
(Input  Signal  Frequency  =  1.0  MHz;  Amplitude  of  the  Signal  =  -30  dB) 
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Figure  9(i).  FFT  Spectral  Response  for  the  DDC  A/D  Convener 
(Input  Signal  Frequency  =  1.25  MHz;  Amplitude  of  the  Signal  =  Full-Scale) 

00 
-10.0 
-20.0 
-30.0 
-40.0 
-50.0 

i  -600 

1  -70.0 

?  -80  0 

2 

-80.0 
-too.o 
-no.o 
-120.0 
-130.0 
-1400 

0.00  0.05  0.10  0.15  0.20  0.25  0.30  0.35  0.40  0.45  0,50 

Inpul  S«gnat  Frequency 

Nomauaad  F,«w.ncy  .  sarr^bnii  Fr«,uancy 

Fi^e  9(j).  FFT  Spectral  Response  for  the  DDC  A/D  Converter 
(Input  Signal  Frequency  =  1.25  MHz;  Amplitude  of  the  Signal  =  -30  dB) 
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Figure  9(k).  FFT  Spectral  Response  for  the  DDC  A/D  Converter 
(Input  Signal  Frequency  =1.5  MHz;  Amplitude  of  the  Signal  =  Full*Scale) 


0.00  0.05  0.10 


Normakjfld  Frequency  - 


0.20  0.25  0.30  0.3S  0.40  0.4$  OJO 

tnput  Signal  Frequency 


Sarrpang  Freqi'ency 


Figure  9(1).  FFT  Spectral  Response  for  the  DDC  A/D  Converter 
(Input  Signal  Frequency  =  1.5  MHz;  Amplitude  of  the  Signal  =  -30  dB) 
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Figure  9(ni).  FET  Spectral  Response  for  the  DDC  A/D  Converter 
(Input  Signal  Frequency  =  2.0  MHz;  Amplitude  of  the  Signal  =  Full-Scale) 
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Figure  9(n).  FFT  Spectral  Response  for  the  DDC  A/D  Converter 
(Input  Signal  Frequency  =  2.0  MHz;  Amplitude  of  the  Signal  =  -30  dB) 
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Table  6.  Measured  Dynamic  Performance  of  12>Bit,  4>MHz  DDC  A/D  Converter 


sm 


Signal 

Frequency 

(MHz) 

Amplitude 

(dB) 

With 

Harmonics 

(dB) 

Without 

Harmonics 

(dB) 

Effective 

Bits 

SFDR 

(dB) 

0.10 

0 

60.6 

63.5 

9.8 

64.7 

0.10 

-1 

60.4 

62.8 

9.9 

65.3 

0.10 

-3 

60.1 

62.1 

10.2 

65.9 

0.10 

-6 

58.0 

59.2 

10.3 

64.5 

0.10 

-12 

53.7 

55.2 

10.6 

63.0 

0.10 

-18 

47.9 

49.5 

10.7 

57.4 

0.10 

-24 

42.1 

43.5 

10.7 

52.6 

0.10 

-30 

35.6 

37.9 

10.6' 

44.2 

0.27 

0 

61.7 

65.9 

10.0 

67.0 

0.27 

-1 

61.5 

65.1 

10.1 

65.4 

0.27 

-3 

60.3 

63.1 

10.2 

65.2 

0.21 

-6 

58.1 

60.3 

10.4 

64.6 

0.27 

-12 

52.7 

54.0 

10.5 

63.3 

0.27 

-18 

47.5 

48.5 

10.6 

58.4 

0.27 

-24 

41.3 

42.2 

10.6 

54.0 

0.27 

-30 

35.0 

36.9 

10.5 

44.1 

0.54 

0 

64.7 

66.1 

10.5 

72.2 

0.54 

-3 

62.6 

63.9 

10.6 

71.6 

0.54 

-6 

60.3 

61.4 

10.7 

70.4 

0.54 

-12 

54.8 

55.6 

10.8 

68.0 

0.54 

-18 

48.5 

49.2 

10.7 

59.0 

0.54 

-24 

42.4 

43.4 

10.7 

53.9 

0.54 

-30 

34.5 

37.1 

10.4 

40.3 

1.00 

0 

62.9 

63.6 

10.1 

76.0 

1.00 

-1 

62.7 

63.5 

10.3 

73.7 

1.00 

-3 

61.0 

61.9 

10.3 

73.4 

1.00 

-6 

58.6 

59.2 

10.4 

72.2 

1.00 

-12 

53.1 

53.5 

10.5 

66.3 

1.00 

-18 

47.6 

48.2 

10.6 

59.2 

1.00 

-24 

41.4 

42.3 

10.6 

51.6 

1.00 

-30 

34.8 

36.7 

10.5 

41.8 
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Table  6.  (Concluded) 


SNR 

Signal 

Frequency 

(MHz) 

Amplitude 

(dB) 

With 

Hannonics 

(dB) 

Without 

Hannonics 

(dB) 

Effective 

Bits 

SFDR 

(dB) 

1.25 

0 

63.1 

65.2 

10.2 

70.8 

1.25 

-1 

62.7 

64.3 

10.3 

72.2 

1.25 

-3 

61.7 

63.3 

10.5 

69.3 

1.25 

-6 

59.7 

61.0 

10.6 

69.6 

1.25 

-12 

54.1 

55.1 

10.7 

62.5 

1.25 

-18 

48.3 

49.2 

10.7 

60.0 

1.25 

-24 

42.3 

43.5 

10.7 

52.6 

1.25 

-30 

35.9 

37.8 

10.7 

45.2 

1.50 

0 

60.4 

62.6 

9.7 

67.6 

1.50 

-1 

61.1 

62.1 

10.0 

70.4 

1.50 

-3 

60.8 

62.4 

10.3 

71.7 

1.50 

-6 

58.7 

60.1 

10.5 

69.4 

1.50 

-12 

53.9 

55.3 

10.7 

64.0 

1.50 

-18 

48.4 

49.5 

10.7 

58.7 

1.50 

-24 

42.4 

43.7 

10.7 

53.9 

1.50 

-30 

36.2 

37.7 

10.7 

45.5 

2.00 

0 

60.1 

60.5 

9.7 

75.6 

2.00 

-1 

59.7 

60.5 

9.8 

71.4 

2.00 

-3 

58.5 

59.7 

9.9 

68.2 

2.00 

-6 

56.9 

57.7 

10.1 

68.2 

2.00 

-12 

52.1 

52.6 

10.4 

64.6 

2.00 

-18 

46.7 

47.0 

10.5 

63.3 

2.00 

-24 

41.0 

41.4 

10.5 

55.9 

2.00 

-30 

35.1 

35.5 

10.5 

50.2 
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signal  amplitudes  near  the  maximum  the  SINAD  is  approximately  equal  to  65  dB  at 
0.54  MHz,  but  at  other  frequencies  it  is  up  to  approximately  6  dB  less. 

Figure  7  shows  that  for  most  signal  amplitudes  and  frequencies,  the  DDC  converter  loses 
approximately  1.3  to  1.5  bits  with  respect  to  an  ideal  12-bit  converter,  which  corresponds  to 
10.5  to  10.7  effective  bits  under  dynamic  conditions.  However,  for  full-scale  amplitude 
signals,  it  loses  approximately  2.2  bits  at  both  the  high  (>1.5  MHz)  and  low  frequency 
(0.1  MHz). 

Hie  SFDR  is  71.0  to  73.0  dB  for  large  signal  amplitudes.  At  other  frequencies  it  is  similar  to 
the  results  at  0.54  and  1.25  MHz,  except  at  very  low  frequencies  of  100  kHz  and  270  kHz, 
where  the  distortion  is  higher  and  the  SFDR  is  65  dB  for  large  signal  amplitudes.  The  B2 
specification  specifies  that  for  a  signal  frequency  that  is  0.54  MHz  near  full-scale,  the 
minimum  RMS  SNR  of  the  converter  must  be  63  dB.  The  SFDR  (for  signal  frequencies  0.54 
and  1.25  MHz)  should  be  70  dB  below  full-scale.  Our  test  results  clearly  show  that  the  DDC 
converter  meets  the  B2  SINAD  and  SFDR  specifications. 

Hgures  10(a)  through  10(g)  show  the  absolute  values  of  the  errors  for  seven  different 
frequencies.  As  shown,  the  error  performance  is  exceptionally  good  at  a  nominal  frequency 
of  544  kHz;  integral  nonlinearity  is  under  +  1.0  LSB.  However,  at  two  lower  frequencies  of 
101  and  271  kHz,  the  third-order  harmonic  distortion  is  noticeable  as  evidenced  by  the 
S-shape  of  the  curve,  although  it  is  higher  and  it  is  within  the  specification.  If  the  noise  is 
visually  disregarded  (which  causes  broadening  of  the  curve),  the  mean  output  error  is  less 
than  approximately  +  1.5  LSB. 


3.2  TEST  RESULTS  FOR  THE  ADC  12.BIT,  4.MHZ  A/D  CONVERTER  -  UNIT  1 

Our  analysis  of  the  test  results  for  the  ADC  12-bit,  4-MHz  A/D  converter  (Unit  1)  (at  two 
different  frequencies)  is  presented  on  the  following  pages,  in  the  form  of  plots,  as  follows: 
noise  versus  input  signal  level  under  dynamic  conditions  (figure  11);  SINAD  versus  input 
*  gnal  level  (figure  12);  bits  lost  versus  input  signal  level  (figure  13);  SFDR  versus  input 
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Figure  10(a).  Errors  of  DDC  A/D  Converter  (Input  Signal  Frequency  =  101  kHz) 


Figure  10(b).  Errors  of  DDC  A/D  Converter  (Input  Signal  Frequency  =  271  kHz) 
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Figure  10(d).  Errors  of  DDC  A/D  Converter  (Input  Signal  Frequency  =  997  kHz) 
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Figure  10(g).  Errors  of  DDC  A/D  Converter  (Input  Signal  Frequency  =  1991  kHz) 
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Figure  1 1 .  Noise  vs.  Input  Signal  Level  for  the  ADC  A/D  Converter  -  Unit  1 


Figure  12.  SINAD  vs.  Input  Signal  Level  for  the  ADC  A/D  Converter  -  Unit  1 
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Figure  14.  SFDR  vs.  Input  Signal  Level  for  the  ADC  A/D  Converter  -  Unit  I 
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signal  level  (figure  14);  and  FFT  spectral  response  for  a  full  amplitude  signal  (figures  15a 
and  ISb). 

Table  7  summarizes  the  measured  dynamic  performance.  The  converter  has  1.15  to 
1.4  LSB  of  RMS  noise  for  full-scale  input  signals.  Compared  to  the  DDC  converter,  it  has 
an  excess  of  0.6  LSB  of  noise.  Figure  13  shows  that  it  loses  2.4  to  2.9  bits  for  full-scale 
input  signals.  However,  it  loses  2.6  to  3.0  bits  for  signals  that  are  30  dB  below  full-scale.  It 
provides  nine  effective  bits  under  dynamic  conditions.  The  SFDR  of  this  converter  is  60  to 
65  dB  for  large  signals. 

The  spectrum  plots  (figures  15a  and  15b)  show  an  increase  of  6  dB  in  noise  floor  over  a 
wideband  in  the  upper  half  of  the  Nyquist  band.  This  may  be  caused  by  high-order 
nonlinearity  which  results  fi’om  the  second-stage  error  correction  circuits  not  performing 
properly. 

Based  on  the  test  results,  the  baseline  ADC  A/D  converter  is  noisier  than  the  DDC  unit  and 
does  not  meet  the  SNR  and  SFDR  specifications.  Since  this  unit  had  such  poor  performance 
at  the  B2  specified  frequencies  of  0.54  MHz  and  1.25  MHz,  respectively,  we  did  not  test  the 
converter  at  five  other  frequencies.  However,  the  baseline  ADC  unit  did  meet  all  the  B2 
requirements  during  qualification  testing  of  the  radar.  Since  it  has  been  in  operation  for  three 
years,  perhaps  its  poor  dynamic  performance  may  be  due  to  some  component  failure  in  the 
board.  In  order  to  come  to  a  final  conclusion,  we  obtained  another  ADC  operational  unit 
from  Paramax  and  then  repeated  the  measurements.  The  next  section  presents  the  details  of 
the  test  results  on  the  ADC  A/D  converter  Unit  2. 


33  TEST  RESULTS  FOR  THE  ADC  12-BlT,  4-MHZ  A/D  CONVERTER  -  UNIT  2 

Our  analysis  of  the  test  results  for  the  ADC  12-bit,  4-MHz  A/D  converter  (Unit  2)  (at  four 
different  frequencies)  is  presented  on  the  following  pages,  in  the  form  of  plots,  as  follows; 
noise  versus  input  signal  level  under  dynamic  conditions  (figure  16);  SINAD  versus  input 
signal  level  (figure  17);  bits  lost  versus  input  signal  level  (figure  18);  SFDR  versus  input 
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Figure  15(a).  FFT  Spectral  Response  for  the  ADC  A/D  Converter  -  Unit  1 
(Input  Signal  Frequency  =  0.54  MHz;  Amplitude  of  the  Signal  =  Full-Scale) 
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Figure  15(b).  FFT  Spectral  Response  for  the  ADC  A/D  Convener  -  Unit  1 
(Input  Signal  Frequency  =  1.25  MHz;  Amplitude  of  the  Signal  =  Full-Scale) 


Table  7.  Measured  Dynamic  Performance  of  12-Bit,  4-MHz  ADC  A/D 

Converter  -  Unit  1 


sm 

Signal  With  Without  Effective 

Frequency  Amplitude  Harmonics  Harmonics  Bits  SFDR 

(MHz)  (dB)  (dB)  (dB)  (dB) 


0.54 

0 

59.4 

61.7 

9.6 

64.4 

0.54 

-1 

58.5 

60.3 

9.6 

64.2 

0.54 

-3 

57.6 

58.9 

9.8 

67.2 

0.54 

-6 

59.9 

56.1 

9.7 

58.6 

0.54 

-12 

46.8 

50.0 

9.5 

50.1 

0.54 

-18 

40.5 

43.2 

9.4 

46.5 

0.54 

-24 

36.2 

39.1 

9.7 

41.5 

0.54 

-30 

28.6 

32.3 

9.4 

32.8 

1.25 

0 

56.5 

60.1 

9.1 

60.7 

1.25 

-1 

56.8 

59.8 

9.3 

62.0 

1.25 

-3 

56.8 

58.1 

9.6 

65.8 

1.25 

-6 

54.2 

56.2 

9.7 

62.0 

1.25 

-12 

49.1 

49.7 

9.9 

67.9 

1.25 

•18 

44.0 

45.2 

10.0 

54.8 

1.25 

-24 

36.1 

39.6 

9.7 

39.3 

1.25 

-30 

25.6 

33.9 

9.0 

26.7 
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Figure  1 8.  Bits  Lost  vs.  Input  Signal  Level  for  the  ADC  A/D  Converter  -  Unit  2 
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Figure  19.  SFDR  vs.  Input  Signal  Level  for  the  ADC  A/D  Convener  -  Unit  2 


signal  level  (figure  19);  and  FFT  spectral  response  for  a  full  amplitude  signal  that  is  30.0  dB 
below  full-scale  at  each  frequency  (figures  20(a)  through  20(h)). 

Table  8  summarizes  the  measured  dynamic  performance.  In  general  the  converter  has 
1.1  LSB  of  RMS  noise  at  most  of  the  frequencies  and  for  different  signal  amplitudes  at  each 
frequency.  However,  at  a  signal  frequency  of  0.1  MHz,  the  noise  rises  appreciably  to 
3.0  LSB,  which  may  be  associated  with  the  missing  code  problem  that  it  experiences  under 
these  conditions.  It  loses  between  2.2  and  2.5  bits  at  all  frequencies  and  different  signal 
amplitudes  at  each  level.  However,  for  a  full-scale  amplitude  signal  with  a  frequency  of 

1.5  MHz,  it  loses  3.5  bits,  which  may  be  again  due  to  the  missing  code  problem.  It  provides 

8.5  to  9.8  effective  bits  under  dynamic  conditions.  The  SFDR  is  62  to  66  dB  for  large 
signals. 

The  spectrum  plots  show  an  increase  of  6  dB  in  noise  floor  over  a  wideband  in  the  upper  half 
of  the  Nyquist  band.  This  may  be  caused  by  high-order  nonlinearity  which  results  from  the 
second-stage  error  conection  circuits  not  performing  properly. 

Figures  21(a)  and  21(b)  show  the  absolute  values  of  the  errors  of  the  converter  for  two 
different  frequencies.  Figure  21(a)  shows  the  converter  has  3.0  LSB  peak-to-peak  of  second 
harmonic  distortion.  The  error  curve  is  not  symmetrical;  it  is  different  for  negative-going 
voltages  (zero  angles  of  zero  to  it)  and  for  positive-going  voltages  (angles  n  to  2;r).  At 
angles  corresponding  to  the  first  sample  after  a  positive  peak  there  are  large  errOTS  up  to 
20  LSB.  On  the  second  conversion  after  the  positive  peak,  the  errors  are  moderately  large 
(approximately  4.0  LSB).  From  the  scatter  of  the  curve  it  can  be  concluded  that  the  RMS 
noise  is  slightly  greater  than  1.0  LSB.  Due  to  the  missing  codes  and  large  errors,  it  has  errors 
greater  than  5.0  LSB  and  does  not  meet  the  integral  nonlinearity  specification. 

Based  on  the  test  results,  even  Unit  2  is  noisier  than  the  DDC  unit  and  does  not  meet  the 
SNR,  SFDR,  and  integral  linearity  B2  specifications. 
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Figure  20(a).  FFT  Spectral  Response  for  the  ADC  A/D  Converter  -  Unit  2 
(Input  Signal  Frequency  =  0.10  MHz;  Amplitude  of  the  Signal  =  Full-Scale) 
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Figure  20(b).  FFT  Spectral  Response  for  the  ADC  A/D  Converter  -  Unit  2 
(Input  Signal  Frequency  =  0.10  MHz;  Amplitude  of  the  Signal  =  -30  dB) 
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Figure  20(c).  FFT  Spectral  Response  for  the  ADC  A/D  Converter  -  Unit  2 
(Input  Signal  Frequency  =  0.27  MHz;  Amplitude  of  the  Signal  =  Full-Scale) 
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Figure  2(Xd).  FFT  Spectral  Response  for  the  ADC  A/D  Converter  -  Unit  2 
(Input  Signal  Frequency  =  0.27  MHz;  Amplitude  of  the  Signal  =  -30  dB) 
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Figure  20(e).  FFT  Spectral  Response  for  the  ADC  A/D  Converter  -  Unit  2 
(Input  Signsd  Frequency  =  0.54  MHj^  Amplitude  of  the  Signal  =  Full-Scale) 
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Figure  20(f).  FFT  Spectral  Response  for  the  ADC  A/D  Converter  -  Unit  2 
(Input  Signal  Frequency  =  0.54  MHz;  Amplitude  of  the  Signal  =  -30  dB) 
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Figure  20(g).  FFT  Spectral  Response  for  the  ADC  A/D  Converter  -  Unit  2 
(Input  Signal  Frequency  =1.5  MHz;  Amplitude  of  the  Signal  =  Full-Scale) 


0.0 

-10.0 

-20.0 

-30  0 

-40.0 

-50.0 

§.  -60,0 

f  -70.0 
■S 

?  -eo.o 

2 

-90.0 

-100.0 

-110.0 

-120.0 

-130.0 

-140.0 


0.00  0.05 


0.10  0l5  0,20 

NormaNzed  Fr«quancy  • 


0.2S  0.30  0.35 

Input  Signal  Frequency 
Saiipiing  Fraquaney 


0  40  0  4S 


0.50 


Figure  20pi).  FFT  Spectral  Response  for  the  ADC  A/D  Converter  -  Unit  2 
(Input  Signal  Frequency  =1.5  MHz;  Amplitude  of  the  Signal  =  -30  dB) 
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Table  8.  Measured  Dynamic  Performance  of  4*MHz  ADC  A/D 

Converter  -  Unit  2 


Signal 

Frequency  Amplitude 
(MHz)  (dB) 


With  Without  Effective 

Harmonics  Harmonics  Bits 

(dB)  (dB) 


0.10 

0.0 

58.8 

60.0 

9.5 

0.10 

-1.0 

57.7 

58.6 

9.5 

0.10 

-3.0 

56.1 

56.8 

9.5 

0.10 

-6.0 

52.2 

52.7 

9.4 

0.10 

-12.0 

46.8 

47.0 

9.5 

0.10 

-18.0 

40.6 

40.9 

9.4 

0.10 

-24.0 

34.8 

35.3 

9.5 

0.10 

-30.0 

28.7 

29.3 

9.5 

0.27 

0.0 

59.2 

60.8 

9.5 

0.27 

-1.0 

58.8 

60.2 

9.6 

0.27 

-3.0 

57.2 

58.6 

9.7 

0.27 

-6.0 

54.4 

55.9 

9.7 

0.27 

-12.0 

48.7 

50.3 

9.8 

0.27 

-18.0 

42.6 

44.0 

9.8 

0.27 

-24.0 

36.2 

38.8 

9.7 

0.27 

-30.0 

29.2 

32.9 

9.5 

0.54 

0.0 

58.1 

60.6 

9.4 

0.54 

-1.0 

58.6 

60.3 

9.6 

0.54 

-3.0 

57.2 

58.4 

9.7 

0.54 

-6.0 

53.5 

55.8 

9.6 

0.54 

-12.0 

46.7 

49.5 

9.5 

0.54 

-18.0 

41.3 

43.7 

9.6 

0.54 

-24.0 

36.1 

37.8 

9.7 

1.50 

0.0 

52.8 

53.6 

8.5 

1.50 

-1.0 

55.2 

56.4 

9.0 

1.50 

-3.0 

55.2 

56.4 

9.4 

1.50 

-6.0 

54.1 

54.9 

9.7 

1.50 

-12.0 

49.3 

50.7 

9.9 

1.50 

-18.0 

44.1 

45.4 

10.0 

1.50 

-24.0 

37.4 

39.0 

9.9 

1.50 

-30.0 

31.4 

32.7 

9.9 

SFDR 

(dB) 


66.5 

66.5 

68.6 
66.6 

65.5 

58.3 

49.6 

43.1 

66.7 

68.3 

67.2 

64.1 

57.8 

54.7 
42.6 
36.0 

62.5 

64.8 

65.1 

60.6 

54.2 
49.0 

44.7 

62.1 

65.3 

66.3 

66.8 

58.5 
54.8 
47.0 

40.6 
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3.4  SUMMARY  AND  CONCLUSIONS  OF  SINE-WAVE  TESTS 


After  comparing  the  dynamic  performance  of  the  ADC  and  DDC  converters,  we  found  that 
the  new  DDC  converter  met  the  integral  linearity,  SINAD,  and  SFDR  B2  specifications  and 
exceeded  the  performance  of  the  two  three- year-old  ADC  units. 
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SECTION  4 


MONOTONICITY  ANALYSIS  AND  TEST  RESULTS 


In  section  3,  we  validated  the  B2  requirements  for  the  RMS-signal-io-RMS-ncise  ratio, 
SFDR,  and  integral  linearity  requirements.  In  this  section  we  validate  the  monoioniciiy 
requirements  of  both  converters.  Based  on  the  spectral  analysis  test  results,  we  concluded 
that  the  two  ADC  converters  were  defective  and  not  functioning  properly.  However,  for 
completeness,  we  present  the  monotonicity  test  results  for  the  ADC  Unit  2  only. 

4.1  TEST  DESCRIPTION 

The  test  is  a  variant  of  the  standard  sine-wave  test.  The  test  signal  is  a  specu-ally  pure  and 
low-noise  filtered  full-scale  sine  wave  at  a  frequency  of  1 25006 1 .03  Hz  { 1 J.  For  the  DDC 
convener,  the  signal  amplitude  is  217 1  LSB  (for  this  convener,  the  LSB  corresponds  to  1 
mv),  sufficient  to  drive  the  A/D  converter  slightly  into  saturation.  The  B2  specification 
requires  measurements  of  all  codes.  This  amplitude  is  the  smallest  that  we  could  generate 
conveniently  which  exercises  all  codes.  For  the  ADC  converter,  the  signal  is  reduced  slightly 
(to  1990  LSB)  to  avoid  saturation  because  the  spectral  analysis  test  indicated  a  possible 
problem  near  saturation,  as  discussed  in  section  3. 

The  sine  wave  is  applied  to  the  A/D  converter  under  test  and  the  resulting  digital  65536 
(64  K)  consecutive  data  points  are  stored  as  a  file  on  a  floppy  disk.  This  is  repeated  as  fast  as 
possible,  with  a  typical  interval  of  about  two  minutes  per  record.  For  each  convener,  a  total 
of  eight  data  records,  representing  512  K  measurements,  are  made.  Final  data  analysis  is 
done  using  a  combination  of  custom  and  commercial  software  running  on  the  mainframe 
VAX  and  PCs. 
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42  DATA  ANALYSIS 


The  first  step  in  the  data  analysis  is  calculating  the  input  cosine  wave  from  the  output  data 
using  a  FORTRAN  orogram,  based  upon  a  sine-wave  curve-fitting  algorithm  described  in 
section  2  [11].  Becj.use  the  convener  is  driven  into  saturation,  the  program  is  modified  to 
eliminate  the  saturated  samples  prior  to  fitting  the  cosine  wave  to  the  rest  of  the  data  (the 
FORTRAN  source  code  is  listed  in  appendix  B).  Each  data  record  of  64  K  sample  points 
represents  16384  microseconds  of  signal.  During  this  interval  the  signal  will  have 
approximately  20481  cycles  (actually  20480.99991552  cycles).  Since  we  have  an  integral 
number  of  cycles  in  the  data  record,  a  perfect  clipped  cosine  wave  is  fitted  to  the  data  with  a 
negligible  error.  For  each  unsaturated  data  point,  the  input  voltage,  its  phase  angle,  and  the 
error  (which  is  the  difference  between  the  estimated  input  and  measured  output)  are 
calculated.  The  data  is  tlien  separated  into  two  groups,  based  upon  the  phase  angle  of  the 
input.  One  group  consists  of  measurements  of  increasing  voltage,  while  the  other  contains 
those  of  decreasing  voltage. 

Figure  10  shows  the  actual  errors  of  the  DDC  converter.  For  most  measurements,  the  mean 
error  is  less  than  +  0.5  LSB,  and  the  variation  from  the  mean  is  also  less  than  +  0.5  LSB. 
Many  subranging  A/D  conveners  behave  differently  with  increasing  or  decreasing  (i.e., 
positive-going  or  negative-going)  signals;  the  software  separates  them  into  distinct  groups 
which  can  be  analyzed  independently,  or  combined  if  desired.  The  two  groups  are  then 
individually  reordered  in  order  of  increasing  input  voltage.  The  input  voltages  and  output 
codes  of  each  group  are  saved  as  disk  files.  A  second  FORTRAN  program  reads  the  eight 
pairs  of  disk  files,  combines  them,  resons  them  in  order  of  increasing  input  voltage,  and  then 
writes  them  to  another  pair  of  data  files.  One  of  these  files  contains  the  input  voltages  and 
output  codes  of  all  256  K  measurements  made  with  increasing  voltage;  the  other  one  contains 
those  made  with  decreasing  voltage. 

Different  records  with  randomly-determined  stan  times,  several  minutes  apan,  will  not  be 
phase-coherent.  The  measured  sample  points  will  not  repeat  exactly  between  records. 
Therefore,  the  combined  eight  records  with  512  K  measurements  are  almost  uniformly 
distributed  in  phase  angle. 
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The  uniform  distribution  of  the  measurements,  with  respect  to  signal  phase  angle,  causes  the 
distribution  of  voltage  measurements  to  vary  inversely  with  the  sine  of  the  phase  angle.  The 
minimum  number  of  measurements  per  code  occurs  at  the  mid-scale  of  the  A/D  converter. 
At  signal-aero  crossing,  there  will  be  approximately  38  measurements^  per  code  made  with 
increasing  voltages,  and  a  similar  number  for  decreasing  voltages;  and  there  will  be  more 
than  100  measurements  per  code  near  the  saturation  points  of  the  converter. 

For  the  case  of  a  noise-free  A/D  converter,  the  code  edges  could  be  determined  to  a 
resolution  of  0.03  (i.e.,  1/38)  LSB  by  using  512  K  data  points,  but  the  resolution  and 
accuracy  is  lower  for  a  real  A/D  converter  with  internal  noise.  To  minimize  the  probability 
of  drawing  erroneous  conclusions  from  the  noisy  data,  we  employ  three  different  methods  to 
evaluate  the  code  widths:  the  histogram,  the  code-mean,  and  the  code-edge  methods. 


4.3  DDC  A/D  CONVERTER  ANALYSIS 
4J.1  Histogram  Method 

The  histogram  is  the  best  possible  test  method  for  detecting  missing  or  wide  codes,  but  it 
does  not  indicate  monotonicity.  If  the  converter  is  known  to  be  monotonic  a  priori,  then  the 
histogram  can  be  used  as  an  indicator  of  the  uniformity  of  the  code  widths.  With  this 
method,  the  derivation  of  the  code  widths  is  based  on  the  number  of  occurrences  of  the  codes. 

A  histogram  of  the  number  of  occurrences  of  each  possible  digital  output  code  is  created. 

The  histogram  is  then  normalized  by  multiplying  by  the  absolute  value  of  the  sine  of  the 
signal  phase  angle  to  compensate  for  the  larger  number  of  observations  expected  near  the 


2  Depends  on  the  number  of  data  samples  and  the  amplitude  of  the  signal.  For  example, 
wiA  512  K  samples  and  the  amplitude  of  the  signal  equal  to  2171.0  LSB, 

the  nuinber  of  measurements  per  code  =  f  |  =  38 

[2l71x2ffJ 
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saturation  points.  After  the  normalization,  the  expected  value  of  the  number  of  occurrences 
of  a  code  is  proportional  to  the  effective  width  of  that  code.  The  normalized  histograms  of 
the  DDC  A/D  converter  are  shown  in  figures  22(a),  22(b),  and  22(c),  for  decreasing  input 
voltage,  increasing  input  voltage,  and  combined  voltages,  respectively.  Horizontal  dashed 
lines  at  the  B2  speciBcadon  limits  of  0.10  and  1.90  dmes  the  nominal  code  width  are  also 
shown,  which  correspond  to  a  minimum  of  four  hits  and  a  maximum  of  72  hits  per  code  in 
our  histogram  analysis. 

The  histogram  shows  that  there  are  no  missing  codes.  For  negative  signal  slope,  there  are  no 
narrow  codes,  and  only  three  codes  are  slightly  wider  than  the  specification  limit  However, 
for  increasing  voltages,  there  is  one  narrow  code  and  seven  wide  codes.  Generally,  most  of 
the  codes  which  are  extremely  wide  or  narrow,  with  signals  slewing  in  one  direction,  behave 
mote  nominally  with  slopes  in  the  opposite  direction.  If  the  two  slope  polarities  are 
combined,  all  code  widths  are  in  the  range  of  0.2  to  2.0  LSB,  and  only  one  code  is  slightly 
wider  than  the  specification. 

4.3.1.1  Sensitivity  Analysis 

The  accuracy  of  the  histogram  method  of  analysis  is  dependent  on  the  number  of  samples  and 
the  noise  of  the  converter,  as  shown  in  appendix  C.  Measurements  made  on  the  DDC 
converter  for  the  spectral  analysis  test  indicate  that  it  has  approximately  0.5  LSB  of  RMS 
noise. 

The  number  of  occurrences  of  each  code  is  affected  by  the  internal  noise  of  the  converter.  In 
appendix  C,  each  number  of  occurrences  is  modeled  as  a  random  variable  with  a  mean  value 
equal  to  the  product  of  the  code  width  (expressed  in  mVoIts  or  LSB)  times  the  density  of  the 
sampling  (expressed  in  samples  per  mVolt).  For  codes  which  are  narrow  compared  to  the 
RMS  noise  level  of  the  converter,  the  number  of  occurrences  approaches  a  Poisson 
distribution,  with  a  variance  equal  to  the  mean.  However,  for  codes  which  are  wide 
compared  to  the  RMS  noise  of  the  converter,  the  variance  in  the  number  of  occurrences 
approaches  a  constant,  which  is  1.13  times  the  product  of  the  RMS  noise  (expressed  in 
mVolts  or  LSB)  times  the  density  of  sampling,  as  shown  in  figure  C-3. 
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Figure  22(b).  Normalized  Histogram  for  the  DDC  A/D  Converter  (Increasing  Voltage) 
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Figure  22(c).  Normalized  Histogram  for  the  DDC  A/D  Converter  (Combined  Slopes) 
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For  the  DDC  converter  the  histogram  for  combined  slopes  (see  figure  22c)  shows  that  the 
sampling  density  is  76  occurrences  per  LSB,  and  that  many  of  the  narrow  codes  have  about 
25  occurrences  corresponding  to  25/76=0.33  LSB,  which  is  less  than  the  0.5  LSB  noise  level 
of  the  converter.  From  figure  C-3,  we  see  that  the  variance  of  the  number  of  occurrences  of 
narrow  codes  is  21  (83  percent  of  the  mean).  The  standard  deviation  of  the  narrow  codes  is 
4.6  observations  (or  4.6/76=0.06  LSB);  hence,  the  measurement  errors  would  be  ±  0. 12  LSB. 

Most  of  the  wide  codes  for  the  DDC  converter  have  about  133  observations,  corresponding  to 
a  width  for  wide  codes  of  133/76  =  1.75  LSB.  Since  this  is  3.5  times  the  (0.5  LSB)  RMS 
noise  of  the  converter,  we  assume  a  normal  distribution  with  a  variance  (see  figure  B-3)  of 
43  =  (0.32  *  133),  The  standard  deviation  of  the  wide  codes  is  6.6  observations  (or 
6.6/76=0.09  LSB).  Hence,  we  can  conclude  from  this  analysis  that  the  largest  observation  of 
2.0  LSB  observed  in  the  DDC  code  analysis  could  have  been  caused  by  random  variation. 
Based  on  the  512  K  samples  used  in  the  histogram  sensitivity  analysis,  it  can  be  concluded 
that  all  code  widths  of  the  DDC  converter  are  1  ±  0.75  LSB,  to  a  measurement  accuracy  of 
±0.15  LSB. 

4.3.2  Code>Mean  Method  of  Analysis 

Proving  that  the  output  of  the  converter  is  a  monotonic  function  of  the  input,  as  defined 
earlier,  would  require  several  hundred  measurements  at  every  possible  input  voltage.  To 
simplify  the  test,  we  employ  a  less  restrictive  definition  of  monotonicity.  For  this  definition 
we  show  that  the  input  is  a  monotonic  function  of  the  discrete  output  code.  For  each  of  the 
observed  output  codes,  we  determine  the  mean  input  voltage  which  produced  it.  If  these 
code-mean  input  voltages  form  a  monotonic  sequence,  then  the  converter  is  considered  to  be 
monotonic.  The  separation  of  the  code-mean  input  voltages  is  then  used  to  estimate  the  code 
widths.  If  the  separation  between  successive  code  means  is  always  positive,  the  converter  is 
monotonic. 

For  this  method  of  the  analysis  we  calculate  the  average  input  voltage  associated  with  each 
possible  output  code.  All  of  the  calculated  code-mean  input  voltages  are  a  monotonic 
function  of  output  code  for  either  slope  polarity.  The  separations  between  successive  values 
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of  the  code  means  are  shown  in  figures  23(a),  23(b),  and  23(c),  for  decreasing  input  voltage, 
increasing  input  voltage,  and  combined  slopes,  respectively.  Horizontal  dashed  lines  at  the 
B2  specification  limits  of  0.10  and  1.90  times  the  nominal  code  width  are  also  shown. 

For  decreasing  voltages,  all  of  the  code  widths  of  the  DDC  converter  are  within  the 
speciHcation,  except  for  one  narrow  code.  However,  for  increasing  voltages,  or  for  the 
combined  case,  all  of  the  code  widths  are  within  specification.  The  minimum  spacing 
between  calculated  code  centers  is  0.05  LSB,  and  the  maximum  spacing  is  1.87  LSB.  When 
both  slope  polarities  are  combined,  the  calculated  linearity  improves.  The  distances  between 
code  centers  are  in  the  range  of  0.1 1  to  1.77  LSB.  Based  on  this  analysis,  we  observe  that  the 
DDC  A/D  converter  is  strictly  monotonic. 

4.33  Code-Edge  Method  of  Analysis 

With  this  method,  we  estimate  the  actual  code  edges  and  their  spacing.  The  lower  edge  of  a 
code  is  defined  as  the  lowest  input  voltage  which  has  at  least  50  percent  probability  of 
producing  that  code  or  higher  at  the  A/D  converter  output  Since  each  of  the  512  K 
measurements  is  taken  at  a  unique  input  voltage,  calculating  the  probability  distribution  of 
the  output  code  is  complicated.  In  our  analysis  we  establish  a  small  window  of  phase  angles 
centered  on  the  code  edge.  For  each  code,  the  window  center  is  adjusted  to  the  lowest  center 
so  that  more  than  half  of  the  occurrences  are  of  that  code  or  higher.  The  width  of  the  sliding 
window  is  arbitrary;  after  evaluating  several  values,  our  analysis  reveals  that  a  window  width 
of  25  gives  optimum  results.  With  the  window  centered  at  the  code  edge,  there  are  12 
measurements  with  an  input  slightly  larger  than  the  edge,  one  at  the  code  edge,  and  12 
slightly  below  it. 

The  separation  between  successive  values  of  the  code  edges  is  shown  in  figures  24(a)  and 
24(b)  for  decreasing  input  voltage  and  increasing  input  voltage,  respectively.  Horizontal 
dashed  lines  at  the  specification  limits  of  0.10  and  1.90  times  the  nominal  code  width  are  also 
shown.  The  measured  width  of  the  wide  codes  is  independent  of  the  choice  of  window 
width.  For  the  positive  slope,  there  are  37  codes  wider  than  1.9  LSB  —  the  widest  is 
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Figure  23(a).  Code-Mean  Separation  for  the  DDC  A/D  Clonverter  (Decreasing  Voltage) 
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Figure  23(b).  Code-Mean  Separation  for  the  DDC  A/D  Converter  (Increasing  Voltage) 
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Figure  23(c).  Code-Mean  Separation  for  the  DDC  A/D  Convener  (Combined  Voltages) 
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Figure  24(b).  Code-Edge  Separation  for  the  DDC  A/D  Converter  (Increasing  Voltage) 
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3.01  LSB.  For  the  negative  slope,  there  are  27  codes  wider  than  1.9  LSB  —  the  widest  is 
2.58  LSB.  Since  there  are  fewer  observations  of  the  narrow  codes,  the  effects  of  noise  are 
greater,  and  the  certainty  of  conclusions  is  lower  and  more  dependent  on  the  choice  of 
window  width.  For  the  positive  slope,  there  are  approximately  six  codes  narrower  than 
0.1  LSB  —  the  narrowest  is  0.007  LSB.  For  the  negative  slope,  there  are  no  codes  that  are 
narrower  than  0.1  LSB;  the  narrowest  is  0.14  LSB.  With  both  slopes  combined,  all  codes  are 
wider  than  0.04  LSB.  Of  the  three  methods  described  for  estimating  code  width,  the  code¬ 
edge  method  is  the  least  accurate,  and  the  most  susceptible  to  noise. 


4.4  ADC  A/D  CONVERTER  MONOTONICITY  ANALYSIS  -  UNIT  2 

The  spectral  analysis  tests  of  the  ADC  Unit  2  indicate  that  it  has  a  serious  problem  when  the 
input  signal  is  a  large  amplitude  sine  wave  of  more  than  1.0  MHz.  Whenever  the  input 
voltage  exceeds  75  percent  of  positive  full-scale,  it  causes  the  next  conversions  to  exhibit 
large  blocks  of  missing  codes  for  negative  voltages  spaced  '28  LSB  apart.  At  1.25  MHz 
there  are  five  blocks  of  more  than  15  consecutive  missing  codes,  and  peak  errors  of  20  LSB. 
Figure  21(a)  shows  the  errors  as  a  function  of  the  input  signal  phase  angle  in  radians  after  the 
positive  peak.  Hence,  prior  to  performing  the  monotonicity  tests,  we  already  know  from  the 
phase-plane  error  analysis  that  the  converter  is  not  strictly  monotonic  and  that  it  has  gross 
errors.  Because  we  suspect  that  the  converter  errors  might  be  related  to  signal  amplitude 
exceeding  positive  saturation,  we  reduce  the  test  signal  amplitude  slightly  (to  1990  LSB), 
which  is  9  LSB  less  than  saturation  because  of  the  converter  offset.  However,  we  use  the 
same  signal  frequency  as  that  used  for  testing  the  DDC  converter.  We  now  present  an 
analysis  using  the  histogram  and  code-mean  methods. 

4.4.1  Histogram  Method 

Figures  25(a),  25(b),  and  25(c)  show  the  normalized  histogram  for  decreasing  voltage, 
increasing  voltage,  and  the  combined  slope  polarities,  respectively.  The  increasing  voltage 
data  has  a  block  of  three  consecutive  missing  codes  which  occur  on  the  third  sample  after  the 
positive  peak.  Figure  25(a)  (decreasing  voltage)  shows  that  there  are  174  missing  codes, 
grouped  into  blocks  that  are  spaced  128  LSB  apart  The  largest  block  is  17  consecutive 
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Figure  25(a).  Normalized  Histogram  for  the  ADC  A/D  Converter  -  Unit  2 
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Figure  25(b).  Normalized  Histogram  for  the  ADC  A/D  Converter  -  Unit  2 
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missing  codes,  and  the  resultant  width  of  the  next  code  is  18  LSB.  By  combining  the 
positive-going  and  negative-going  data,  no  missing  codes  are  found. 


4.4.2  Code-Mean  Method 

Figures  26(a),  26(b),  and  26(c)  show  the  code  widths  computed  as  the  difference  between 
successive  code  means  for  decreasing  voltage,  increasing  voltage,  and  the  combined  slope 
polarities,  respectively.  For  decreasing  voltage,  near  the  edge  of  the  missing  codes  some 
codes  are  not  completely  missing,  but  have  an  abnormally  low  number  of  observations.  The 
code-mean  voltages  calculated  for  these  codes  form  a  non-monotonic  sequence,  with  a 
minimum  code  separation  of  -1.15  LSB.  By  combining  the  positive-going  and  negative- 
going  data,  no  missing  codes  are  found.  The  mean-code  voltages  form  a  non  monotonic 
sequence,  with  a  minimum  spacing  of -8.90  LSB,  and  a  maximum  spacing  of  10.64  LSB. 

This  converter  illustrates  the  importance  of  processing  the  data  according  to  the  phase  of  the 
input  signal,  thereby  separating  the  positive-going  measurements  from  the  negative-going 
ones.  Had  we  combined  both  slope  polarities,  there  would  not  have  been  any  missing  codes, 
and  we  would  not  have  observed  an  obvious  explanation  for  the  exuemely  nonlinear 
behavior. 


4.5  SUMMARY  AND  CONCLUSIONS  OF  MONOTONICITY  ANALYSIS 

Table  9  summarizes  the  monotonicity  analysis  of  the  converters  using  the  three  different 
methods  (the  histogram,  code-mean,  and  code-edge  methods).  Since  many  subranging  A/D 
converters  behave  differently  with  increasing  or  decreasing  (i.e.,  positive-going  or  negative¬ 
going)  signals,  we  present  the  analysis  for  positive-going  or  negative-going,  and  combined 
slopes  to  gain  better  insight  Based  on  our  analysis,  we  conclude  that  the  DDC  converter  is 
strictly  monoionic  and  marginally  meets  the  UAR  B2  specification  requirements.  As 
indicated  in  table  2,  the  DDC  converter  appears  to  have  several  codes  which  may  be  slightly 
wider  than  the  B2  specification.  Confirming  this  with  high  confidence  would  require  many 
more  samples  than  we  used  in  our  analysis  (greater  than  512  K).  The  ADC  converter  is 
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Figure  26(b).  Mean  Code  Width  for  the  ADC  A/D  Converter  -  Unit  2 

(Increasing  Voltage) 
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Figure  26(a).  Mean  Code  Width  for  the  ADC  A/D  Converter  -  Unit  2 
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Figure  26(c).  Mean  Code  Width  for  the  ADC  A/D  Convener  -  Unit  2 

(Combined  Slopes) 
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nonmonotonic  and  has  gross  enors,  particularly  for  negative-going  voltages,  and  it  does  not 
meet  UAR  B2  specification  requirements. 


Table  9.  Monotonidty  •  Diflferentiai  Nonlinearity  Test  Results 


DDC  A/D  Converter 


Test  Method 

Converter 

Narrower 
Than  Spec 

Narrowest 

Code 

iaILSB) 

Widest 

Code 

io.,(.L5B) 

Specification 

DDC 

0 

0.10 

1.90 

Histogram  E)ecieasing 

DDC 

0 

0.22 

2.38 

Histogram  Increasing 

DDC 

I 

0.05 

2.51 

Histogram  Combined 

DDC 

0 

0.31 

2.01 

Code  Mean  Decreasing 

DDC 

1 

0.05 

1.87 

Code  Mean  Increasing 

DDC 

0 

0.11 

1.77 

Code  Mean  Combined 

DDC 

1 

0.04 

1.75 

Code  Edge  Decreasing 

DDC 

0 

0.14 

2.58 

Code  Edge  Increasing 

DDC 

6 

0.007 

3.00 

Best  Estimate 

DDC 

0 

0.2510.15 

1.7510.15 

Test  Method 

Converter 

IXI2£ 

ADC  A/D  Converter 

Narrowest 

Narrower  Code 

IhanSpCg  in(LSB) 

Widest 

Code 

in  (LSB) 

Specification 

ADC 

0 

0.10 

1.90 

Histogram  Decreasing 

ADC 

174 

0.00 

20.92 

Histogram  Increasing 

ADC 

3 

0.00 

2.64 

Histogram  Combined 

ADC 

0 

0.30 

11.22 

Code  Mean  Decreasing 

ADC 

174 

-1.15 

10.11 

Code  Mean  Increasing 

ADC 

3 

-2.50 

2.70 

Code  Mean  Combined 

ADC 

0 

-8.90 

10.64 

Wider 
Than  Spec 

0 


76 


fTi  I—— .000t''r''0 
<N  m 


SECTION  5 


TRANSIENT  AND  OVERVOLTAGE  RECOVERY  TESTING  RESULTS 


In  sections  3  and  4,  we  validated  the  B2  requirements  for  the  RMS-signal-to-RMS  noise 
ratio,  SFDR,  and  monotonicity  requirements.  The  B2  specification  requires  that  the  A/D 
converter  recover  from  any  transient  or  overrange  input  condition  to  an  accuracy  of 
+  5.0  mV  within  200  nanoseconds.  The  converter  should  be  tested  for  all  baseline  voltages 
within  its  ±  2.047  V  linear  range,  with  pulses  whose  peak  voltages  span  the  allowable 
+  4.0  V  overvoltage  input  range.  In  this  section,  we  validate  the  transient  and  overvoltage 
recovery  requirements  of  both  converters. 


5.1  TEST  DESCRIPTION 

The  test  is  a  variant  of  the  standard  sine-wave  test  described  in  section  2.  The  input  signal 
applied  to  the  A/D  converter  under  test  is  the  sum  of  three  signal  components:  a  (dc)  level 
which  is  close  to  the  saturation  point  of  the  converter;  a  small  70-Hz  triangular  modulation 
added  to  the  dc  level;  and  a  train  of  large  amplitude  pulses  superimposed  on  the  two  signal 
components.  The  pulses  are  one  microsecond  wide,  repeated  with  a  period  of  four 
microseconds.  The  input  signal  is  applied  to  the  A/D  converter  under  test  and  the  resulting 
digital  output  is  stored  in  the  Tektronix  High-Speed  Data  Cache  (a  large,  fast,  sequentially- 
addressed  memory),  and  stored  as  a  file  on  a  floppy  disk. 

Since  the  pulse  repetition  period  is  16  cycles  of  the  4-MHz  clock,  a  data  record  of  64  K 
measurements  includes  4096  pulses.  In  the  absence  of  external  triangular  modulation,  the 
measurements  repeat  exactly,  except  for  the  converter  and  test  fixture  noise.  With  a  one- 
microsecond  nominal  pulsewidth,  there  are  four  measurements  of  the  pulse,  followed  by 
twelve  measurements  of  the  settling  between  pulses.  By  applying  a  small  70-Hz  triangular 
modulation,  a  full  cycle  of  the  modulation  is  obtained  in  one  data  record,  allowing 
investigation  of  the  converter  performance  over  a  limited  range  of  input  voltages.  Note  that 
in  this  rest  the  voltages  of  each  pulse  are  modulated,  but  the  relative  phase  of  each 
measurement  in  each  pulse  is  the  same. 
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5.2  TEST  CIRCUIT  DESCRIPTION 


As  discussed  in  section  1,  we  designed  and  built  a  test  fixture  which  can  generate  a  pulse  to 
12  bits  accuracy.  The  schematic  diagram  of  the  fixture  is  shown  in  figure  27,  and  the 
assembly  print  is  shown  in  figure  28.  The  details  of  the  test  fixture  arc  discussed  below. 

In  normal  operation,  a  square  transistor-transistor-logic  (Til,)  level,  4-MHz  clock  is  applied 
to  the  test  circuit  The  input  clock  is  buffered  by  integrated  circuit  2L8,  and  provided  as  an 
output  signal  to  drive  the  A/D  converter.  A  counter  circuit  (25  and  Z6)  driven  by  the  input 
clock  produces  a  digital  pulse  that  is  one  microsecond  (four  clock  cycles)  wide,  with  a 
frequency  of  250  kHz.  The  digital  pulse  is  delayed  by  latch  Z7  for  an  interval  which  is 
adjustable  via  potentiometer  P2,  then  applied  to  a  current  switch  formed  by  transistors  Q1  to 

The  current  switch  superimposes  an  analog  voltage  pulse  upon  a  baseline  voltage  which  is 
established  by  operational  amplifier  Zl.  Switches  SWl  and  SW2  control  the  polarity  and 
amplitude  of  the  analog  pulse,  respectively.  The  baseline  voltage  consists  of  a  dc  component 
which  is  set  via  potentiometer  PI,  plus  an  optional  small  modulation  which  is  an  inverted  and 
attenuated  replica  of  a  low-frequency  signal  applied  to  the  external  modulation  input.  The 
current  switch  is  designed  so  that  the  load  on  Zl  is  independent  of  the  state  of  the  switch. 

The  baseline  voltage  at  the  output  of  amplifier  Zl  is  bypassed  to  system  ground  by  a  large, 
low-impedance  capacitor  network.  This  prevents  it  from  changing  rapidly,  and  insures  that 
the  output  voltage  is  exceptionally  clean  and  stable  in  the  time  interval  between  pulses.  The 
resultant  pulse  rise  and  peak  are  only  moderately  stable,  but  the  fall  to  baseline  is  fast  and 
clean,  with  no  ringing  or  overshoot  The  pulse  fall  time  is  completely  dominated  by  a  single 
(resistor-capacitor  (RC))  time  constant.  The  effective  resistance  is  250  ohms,  and  the  total 
capacitance  is  approximately  50  picoFarads,  which  is  mainly  from  the  short  cable  used  to 
coimect  the  test  fixture  to  the  A/D  converter.  The  output,  therefore,  has  a  time  constant  of 
approximately  13  nanoseconds.  This  means  that  a  4  V  pulse  settles  to  within  ±  5.0  mV  in 
87  nanoseconds,  and  to  ±  2.0  mV  in  99  nanoseconds.  Smaller  pulses  settle  slightly  faster. 
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Figure  27.  Pulse  Generator  Schematic 


5.3  DATA  ANALYSIS 


The  64  K  measured  output  codes  are  entered  into  a  computer  spreadsheet,  and  reordered  to 
fill  a  table  of  4096  lines  of  16  columns.  By  reading  horizontally  along  a  line,  we  get  16 
sequential  measurements  of  one  pulse.  By  reading  the  table  vertically  down,  we  get 
measurements  made  on  4096  pulses,  at  the  same  relative  phase. 

5.3.1  DDC  A/D  Converter  Overvoltage  and  Tranaent  Recovery  Results 

Figures  29  and  30  show  the  performance  of  the  DDC  A/D  converter  with  baseline  voltages 
ranging  from  -1.9  to  -2.05  V  near  the  negative  saturation,  with  a  150  mV  triangular 
waveform  superimposed.  A  +  4.5  V  pulse  drives  the  converter  into  positive  saturation. 

Using  this  modulation,  the  baseline  voltage  changes  by  only  0.08  mV  per  pulse,  so 
effectively  we  are  testing  at  all  input  voltages  over  a  150  mV  range.  In  both  figures  29  and 
30,  the  horizontal  axis  represents  the  number  of  the  pulse  after  the  start  of  the  record,  i.e., 
four  microseconds  per  pulse.  Figures  29  and  30  differ  slightly  (70  nanoseconds)  in  the  delay 
between  the  pulse  and  the  A/D  converter  clock.  In  our  test,  we  take  data  for  two  different 
values  of  circuit  delays  between  the  pulse  and  the  A/D  converter  clock  of  20  and  90 
nanoseconds.  We  choose  those  two  delays  so  that  the  time  available  until  the  next  measured 
sample  of  the  pulse  in  one  case  is  slightly  more  than  200  nanoseconds,  and  in  another  case  it 
is  slightly  less  than  200  nanoseconds.  In  figure  29,  the  nominal  start  of  the  pulse  occurs 
20  nanoseconds  after  a  clock,  and  230  nanoseconds  before  the  next  A/D  converter  clock 
triggers  another  measurement  In  figure  30,  the  voltage  pulse  is  delayed  to  90  nanoseconds 
after  a  clock,  so  it  is  only  160  nanoseconds  before  the  next  A/D  converter  clock  triggers  the 
next  measurement 

Each  figure  has  three  triangular  waves  showing  the  converter  output  codes.  The  highest 
triangle  is  the  first  measurement  after  the  end  of  the  pulse.  The  middle  triangle  is  the  next 
measurement,  and  the  lowest  triangle  is  the  last  measurement,  just  prior  to  the  next  pulse.  In 
both  figures  29  and  30,  the  triangular  outputs  are  smooth  and  have  uniform  separation,  which 
indicates  that  the  converter  has  no  major  nonlinearities  in  this  region,  and  that  the  settling 
response  is  essentially  independent  of  the  final  voltage.  The  bottom  parts  of  figures  29  and 
30  show  the  spread  between  the  largest  and  smallest  values  of  the  1 1  measurements  between 
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Figure  29.  Overvoltage  Recovery  Response  of  DDC  A/D  Converter  When  Pulsed  into 

Positive  Saturation  (20  nsec  Delay) 


Figure  30.  Overvoltage  Recovery  Reponse  of  DDC  A/D  Converter  When  Pulsed  into 

Positive  Saturation  (90  nsec  Delay) 


consecutive  pulses.  In  Hguie  29  it  is  9  ±  1  mV;  however,  in  figure  30  it  is  17  ±  1  mV,  since 
the  pulses  have  less  time  to  settle. 

Figures  31  and  32  show  the  settling  response  on  an  expanded  scale,  when  the  converter  is 
pulsed  into  positive  saturation  for  circuit  delays  of  20  and  90  nanoseconds,  respectively. 

Each  figure  shows  the  output  codes  for  1 1  consecutive  measurements  of  the  settling  interval 
between  pulses.  Three  lines  on  each  figure  show  the  output  at  intervals  separated  by  200 
microseconds  (i.e.,  50  pulses,  or  about  4  mV  in  input  voltage).  The  bottom  curve  in  figure  31 
shows  the  output  code  against  the  compressed  full-scale  range  axis  on  the  right  edge  of  the 
figure.  On  this  scale,  sample  number  zero  is  made  at  an  input  of  about  +  300  LSB  on  the 
falling  edge  of  a  pulse,  and  sample  number  12  is  made  at  about  the  -900  LSB  level  of  the 
next  pulse.  The  other  samples  not  shown  are  all  at  positive  saturadon,  i.e.,  +  2047  LSB. 

With  an  ideal  infinitely  fast  pulse  with  zero  fall  time,  we  would  like  all  of  the  measurements 
between  pulses  to  produce  the  same  output  code.  However,  because  the  fall  dme  is  not 
infinitely  fast  and  there  are  delays  in  the  converter,  the  actual  output  code 
exponentially  to  its  final  value.  For  the  case  where  the  measurement  starts  230  nanoseconds 
after  the  start  of  the  fall  of  the  pulse  (i.e.,  with  a  20-microsecond  delay),  a  sample  is  taken 
near  the  zero  crossing  of  the  signal,  and  sample  12  is  taken  near  the  middle  of  the  rise  of  the 
next  pulse.  This  shows  that  sample  1  is  250  nanoseconds  after  the  midpoint  of  the  falling 
edge  of  the  pulse.  As  discussed  earlier  (from  figure  3),  the  output  has  settled  to  9  mV  in  an 
additional  250  nanoseconds,  whereas  the  B2  specifies  a  settling  accuracy  of  ±  5  mV  within 
200  nanoseconds.  Figure  30,  showing  a  settling  time  of  250  -  70  =  180  nanoseconds  from  the 
midpoint  of  the  falling  edge  of  the  pulses,  has  17  LSB  of  decay. 

It  can  be  concluded  from  our  analysis  that  the  DDC  A/D  converter  does  not  meet  the 
overvoltage  and  transient  recovery  radar  B2  requirements. 

5.3.2  ADC  A/D  Converter  Overvoltage  and  Transient  Recovery  Results  -  Unit  2 

Figures  33  and  34  show  the  performance  of  the  ADC  Unit  2  with  baseline  voltages  ranging 
from  -1.7  to  -2.02  V  near  the  negative  saturation,  with  a  320  mV  triangular  waveform 
superimposed.  There  is  a  4.5  V  amplitude  pulse  driving  the  converter  into  positive 
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Sample  Number  (Time  -  Sample  Number  x  0.25  psec) 

Figure  31,  Overvoltage  Recovery  Reponse  of  DDC  A/D  Converter  When  Pulsed  into 
Positive  Saturation  on  an  Expanded  Scale  (20  nsec  Delay) 
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Figure  32.  Overvoltage  Recovery  Response  of  DDC  A/D  Converter  When  Pulsed  into 
Positive  Saturation  on  an  Expanded  Scale  (90  nsec  Delay) 
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saturation.  Using  this  moduiation,  the  baseline  effectively  changes  by  only  0.16  mV  per 
pulse,  so  we  are  testing  at  all  input  voltages  over  a  320  mV  range.  Figure  33  shows  that  this 
converter  has  a  serious  problem,  since  the  output  is  not  a  clean  triangular  waveform.  At 
certain  voltages,  there  are  large  groups  of  missing  codes  spaced  128  LSB  apart.  The 
behavior  in  one  of  the  missing  code  regions  is  substantially  different  from  that  in  the  good 
regions.  The  peak  errors  are  up  to  40  LSB.  Even  three  microseconds  after  the  voltage  has 
stabilized  (i.e.,  reached  the  unsaturated  voltage  level)  the  output  codes  for  this  converter  are 
in  error  and  certain  codes  do  not  occur. 

Figure  34  shows  an  expanded  view  of  the  settling  for  three  different  input  voltages.  The 
upper  two  curves,  which  are  in  a  good  input  voltage  region,  show  rapid  settling  to  within 
5  LSB  of  the  final  value.  However,  the  lowest  curve,  at  the  edge  of  a  missing  code  region, 
shows  interesting  settling  behavior,  as  indicated  in  figure  34.  The  converter  starts  to  function 
properly  on  samples  1  and  2,  but  on  sample  3  it  suddenly  produces  an  error  of  -30  LSB.  It 
then  recovers,  and  works  correctly  for  samples  4  through  7.  Samples  8  and  9,  which  are  two 
microseconds  after  the  end  of  the  pulse,  again  produce  an  error  of  -30  LSB.  Because  the 
settling  problem  persists  for  more  than  three  microseconds  after  the  end  of  the  pulse, 
changing  the  pulse  delay  by  another  70  nanoseconds,  as  done  for  the  DDC  convener,  will  not 
reveal  any  additional  information.  Therefore,  we  test  this  converter  only  for  a  20- 
nanosecond  delay  between  the  pulse  and  the  A/D  converter  clock.  A  comparison  of  the  DDC 
and  ADC  transient  overvoltage  recovery  characteristics  clearly  indicates  mat  the  DDC 
converter  has  a  better  performance. 


5.4  BEAT  FREQUENCY  TEST 

In  section  5.3,  we  present  the  settling  response  of  the  converters  with  a  time  resolution  of 
0.25  microsecond  per  measurement  By  combining  many  different  data  records  taken  with 
different  circuit  delay  times  set  by  the  potentiometer  P2,  it  would  be  possible  to  determine 
the  settling  response  to  a  better  time  resolution  by  taking  many  measurements  with  various 
circuit  delays  and  combining  the  different  data  records  in  software.  However,  the  beat 
frequency  test  can  be  used  to  easily  determine  the  settling  response  to  a  fine  time  resolution. 
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This  test  method  allows  us  to  obtain  a  detailed  record  of  the  settling  response  of  the  convener 
for  one  particular  input  voltage.  In  the  test  method  described  previously,  we  show  that  the 
settling  response  of  the  DDC  converter  is  independent  of  the  input  voltage.  Hence,  this 
measurement  technique  at  one  voltage  is  representative  of  all  other  voltages  near  the  negative 
saturation.  In  this  section  we  present  beat  frequency  lest  results  for  transient  recovery  and 
overvoltage  recovery  time  of  the  converter. 

The  first  two  cases  are  for  the  transient  recovery  specification.  The  dc  baseline  voltage  is  set 
in  the  linear  range  near  the  saturation  point  (negative  and  positive  voltages,  respectively),  and 
then  it  is  pulsed  with  a  signal  that  drives  it  near  the  opposite  saturation  point  We  also 
perform  additional  tests  which  drive  the  converter  into  saturation  in  order  to  measure  the 
overvoltage  recovery  time.  The  next  section  describes  the  details  of  the  beat  frequency  test 
setup. 


5.4. 1  Test  Description 

The  beat  frequency  test  uses  two  clocks,  one  for  the  pulse  generator  test  equipment,  and 
another  for  the  A/D  converter.  They  are  at  slightly  different  frequencies,  so  their  relative 
phases  change  slowly,  allowing  us  to  measure  all  phases  of  the  pulse  waveform. 

For  this  test,  the  A/D  converter  is  clocked  at  4.0  MHz  by  a  Hewlett  Packard  pulse  generator. 
A  second  pulse  generator  is  phase-locked  to  the  first,  and  set  to  produce  a  slightly  higher 
frequency  (4(X)3910.07  Hz).  This  second  pulse  generator  is  then  used  to  drive  the  test  circuit, 
which  produces  the  signal  voltage  pulses  at  250244.38  pulses  per  second.  This  frequency  is 
chosen  so  that  the  tester  receives  1024  clocks  in  the  time  that  the  A/D  converter  receives 
1023  clocks,  and  the  phase  of  the  signal  being  sampled  will  slowly  increase.  For  example, 
assume  that  at  time  zero,  both  clocks  are  in  phase  and  samples  are  being  taken  every  0.25 
microsecond.  After  3.99609375  microseconds,  the  test  circuit  has  received  16  clocks  and  has 
produced  exactly  one  cycle  of  the  pulse  wavefonns  so  that  the  waveshape  repeats  every 
3.99609375  microseconds.  However,  the  sixteenth  sample  is  not  taken  until 
(4-3.9960375  microseconds)  =  3.40625  nanoseconds  after  we  start  the  second  pulse.  At  this 
instant  the  second  pulse  signal  amplitude  is  equal  to  the  amplitude  of  the  first  pulse  at  3.9 
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nanoseconds.  When  the  tester  receives  1024  clocks,  it  produces  64  pulses,  and  the 
converter  makes  1023  measurements,  so  it  is  once  again  in  phase  with  the  tester  ckKk.  in  the 
first  1023  microseconds  of  a  1024  microsecond  data  record,  which  includes  4092  A/D 
converter  measurements  of  256  pulses,  the  relative  phase  of  the  convener  and  pulse  clocks 
slowly  changes  by  four  clock  periods  (or  one  microsecond).  Reordering  the  sequence  of 
measurements  is  equivalent  to  measuring  four  complete  cycles  of  the  pulse  waveform  with  a 
time  resolution  of  3.90625  nanoseconds. 

5.4J1  DDC  A/D  Converter  Results  from  the  Beat  Frequency  Test 
5.4J1.1  Transient  Recovery  Performance 

Figure  35  shows  the  pulse  response  of  the  DDC  A/D  convener  with  a  large  positive  pulse, 
which  causes  voltage  swings  from  a  baseline  near  the  negative  limit  of  - 1 .75  V  to  near  the 
positive  limit  of  +  1.75  V.  The  output  waveform  is  extremely  clean,  with  no  visible  ringing  or 
overshoot 

Figure  36  is  an  enlarged  display  of  the  initial  region  from  zero  to  300  nanoseconds  of  the 
settling  response  at  the  end  of  the  pulse.  The  output  has  settled  within  5  LSB  of  its  final 
value  in  280  nanoseconds  after  the  end  of  the  pulse.  Equivalently  the  output  has  settled 
within  5  LSB  of  its  final  value  in  240  nanoseconds  after  the  pulse  has  dropped  to  mid-scale. 

If  we  calculate  starting  time  from  when  the  input  pulse  has  reached  5  LSB  of  its  final  value 
(87  nanoseconds),  the  converter  only  requires  an  additional  193  nanoseconds,  which  would 
be  within  specification.  Since  the  dc  baseline  voltage  is  set  in  the  linear  range,  it  can  be 
concluded  that  the  DDC  converter  recovers  from  the  transient  input  condition  in  either  193, 
240,  or  280  nanoseconds  within  5  LSB  of  its  final  value;  however,  the  B2  specification 
requires  that  it  recover  in  200  nanoseconds. 

5  4JL2  Overvoltage  Recovery  Performance 

Figure  37  shows  the  overvoltage  recovery  performance  of  the  DIXT  converter  when 
measuring  large  negative  voltages  and  recovering  from  a  previous  measurement  in  a  positive 
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Figure  35.  Pulse  Response  of  DDC  A/D  Cbnvertcr  with  a  Negative  Base  Voltage 

and  a  1  psec  Pulse  of  -f  3.5V 


Figure  36.  Settling  Response  of  DDC  A/D  Converter  with  a  Negative  Base  Voltage 
and  a  1  ^sec  Pulse  of  4^3.5  V  on  an  Expanded  Scale 
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saturation.  In  this  test  the  converter  recovers  from  the  overrange  input  condition  in 
361  nanoseconds  from  the  first  unsaturated  measurement,  within  5  LSB  of  its  final  value. 

Figure  38  shows  the  overvoltage  recovery  performance  of  the  DDC  convener  when 
measuring  large  negative  voltages  and  recovering  from  a  previous  measurement  in  a  hard 
negative  saturation.  In  this  case,  it  recovers  in  only  155  nanoseconds  after  it  comes  out  of 
saturation,  within  5  LSB  of  its  final  value. 

Figure  39  shows  the  pulse  response  of  this  converter  for  a  signal  with  a  dc  value  set  near  the 
positive  saturation  limit  of  the  converter  and  pulsed  with  a  -3.5  V  negative  amplitude  puLse  to 
near  the  negative  limit  of  the  convener.  Figure  40  shows  the  recovery  response  of  this 
converter  on  an  expanded  scale  recovering  from  a  large  negative  -4.5  V  amplitude  pulse 
which  drove  the  converter  into  negative  saturation.  The  converter  recovers  in  354 
nanoseconds  after  the  end  of  saturation,  within  5  LSB  of  its  final  value.  We  can  model  the 
A/D  converter  overvoltage  recovery  response  as  a  simple  RC  circuit  with  a  time  constant  of 
approximately  55  nanoseconds.  It  can  be  concluded  from  the  analysis  that  the  DDC 
converter  does  not  meet  the  B2  overvoltage  recovery  requirement  of  200  nanoseconds;  a  time 
constant  of  less  than  30  nanoseconds  would  be  required. 


5.5  SUMMARY  AND  CONCLUSIONS  OF  TRANSIENT  AND  OVERVOLTAGE 
RECOVERY  TESTS 

Table  10  compares  overvoltage  and  transient  recovery  characteristics  of  the  DDC  A/D 
converter  and  the  ADC  Unit  2.  The  characteristics  of  the  DDC  converter  are  better  than 
those  of  the  ADC  converter  currently  in  use;  however,  it  does  not  meet  the  transient  and 
overvoltage  NWS  B2  specifications.  In  this  study,  we  did  not  evaluate  the  impact  that  this 
will  have  on  the  overall  system  performance.  We  recommend  that  the  final  step  in  the 
evaluation  be  to  only  reevaluate  the  impact  that  this  will  have  on  the  NWS  performance. 
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Bgure  37.  DDC  A/D  Converter  Overvoltage  Recovery  Response  from  a  +4.8 V  Pulse 


-2020  -  ■  /  - -  -  '  . -  -  ^  . —  - -  ■ 

-2025  -.  !  -  ■  '■ 

-2030  -  j  . - 

-2035  ~  >  . 

-2040 

-2045  -  . .  .  „  .. 

-2050  - , - ^ - 1 - ; - - - 

0  SO  100  150  200  250  300 

Time  (nsec) 


Figure  38.  DDC  A/D  Converter  Overvoltage  Recovery  Response  from  a  -1.7  V  Pulse 
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Figure  39.  DDC  A/D  Converter  Pulse  Response  with  a  Positive  Base  Voltage 

and  a  1  psec  Pulse  of  -3.5V 


I 

§ 


I960 


0  50  too  150  200  250  300  350  400  450  500 

Time  (nsec) 


Figure  40.  DDC  A/D  Converter  Response  When  Recovering  firom  a  -4.4V  Pulse 
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TaUe  10.  Overvoltage  and  Transient  Recovery  Performance  Comparison  of  Two  AJD 

Converters  for  NWS 


A/P  Convener  Typg 


B2 


Parameters 

Specification 

APC 

DDC 

Transient  errors  (in  mV) 

(0.2  microsecond) 

5 

40 

7 

Overvoltage  Recovery  Error 
(in  mV)  (0.2  microsecond) 

5 

50 

12 

Transient  recovery  time 
to  S  LSB  of  its  final  value 
(in  microseconds) 

0.2 

>3 

0.24  -  0.28 

Overvoltage  recovery  time 
to  5  LSB  of  its  final  value 
(in  microseconds) 

0.2 

>3 

0.34 
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APPENDIX  A 

GENERALIZED  SINE-WAVE  ESTIMATION  ALGORITHM 
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A.l  SINE-WAVE  ESTIMATION  ALGORITHM 


Given  samples  of  a  waveform  which  is  dominated  by  a  sine  wave  with  small  additional 

harmonic  distortion  plus  noise,  it  is  desirable  to  develop  an  algorithm  to  estimate  the  signal 
and  noise  parameters  [11].  We  start  with  a  set  of  N  measurements,  V''^  ,  for  i~  1,2 . N  of  a 

waveform  consisting  of  a  large  fundamental  component,  plus  measurement  errors  and 
noise,  e^.  The  fundamental  is  determined  by  four  parameters:  the  offset,  M ;  amplitude,  A ; 

radian  frequency,  tu;  and  initial  phase,  6. 


=  M-¥  A  sin(ty/+  0) 

i  =  l,2..,N 

(A-l) 

=  F-+€- 
1 

i  =  l,2...N 

(A-2) 

where  is  the  error  in  each  measurement  caused  by  noise  and  the  higher  order  harmonics. 
Given  the  V^'s,  we  derive  estimators  for  the  four  parameters  of  the  fundamental,  plus  the  root 
mean  square  (RMS)  value  of  the  measurement  noise,  e,  .  We  establish  an  initial 

approximation  of  the  parameters,  and  then  refine  them  by  performing  several  iterations  of  an 
algorithm  that  minimizes  the  mean  squared  ~rror.  After  iteration  j ,  the  parameter  estimators 
are  Mj,  Aj,  (Oj  and  Oj.  The  speed  and  convergence  of  this  iterative  algorithm  are 

dependent  on  the  initial  estimate. 

A.1.1  Initial  Estimates 

The  initial  estimates  of  the  offset  and  amplitude  arc  obtained  from  the  mean  and  standard 
deviations  of  the  samples. 


1=1 


(A-3) 


(A-4) 
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(A-4) 


I2  ^  ^ 

Once  the  initial  estimate  of  the  offset  is  found,  each  of  the  data  points  is  ex  pined  to 
determine  whether  it  is  above  or  below  the  offset.  The  first  pair  of  data  poii  i  which  are  on 
opposite  sides  of  Mq  determines  a  pair  of  i  values.  iF  and  if +  1.  The  rest  of  uh  ..ata  is 

scanned  to  determine  NC,  the  number  of  complete  signal  cycles  in  the  data  record.  The  last 
offset  crossing  determines  a  pair  of  /  values,  iL  and  iL  + 1 .  The  pairs  of  data  points  can  be 
linearly  interpolated  to  find  a  zero  order  estimate  of  the  time  of  the  first  and  last  offset 
crossings,  FIRSTq  and  LASTq. 


^iF  ~  V+1 


LAST„  =  iL+ 

”  ''iL-ViL^X 


The  initial  estimator  for  the  radian  frequency  is  Q)q. 


(A-5) 


(A-6) 


0}q  =  InNQiLAST  -  FIRST)  (A-7) 

The  parameter  estimates  are  then  used  to  fit  a  segment  of  a  sine  wave  through  the  pair  of 
points  surrounding  the  first  offset  crossing,  and  then  to  solve  for  the  offset  crossing  of  that 
sine  wave,  FIRSTq .  In  a  similar  manner,  an  improved  estimate  of  the  last  offset  crossing, 
LASTp  is  obtained.  Tbese  improved  crossing  estimates  are  then  used  in  equation  (A-7)  to 
obtain  an  improved  estimate  of  the  radian  frequency  G)|. 
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FIRST^  =  if'  -f  WQsin  ^ 


V-Mq 


/ 


L45Tj  = 


(A-8) 


where 


iF'  = 


iF 

< 

/or|l^f-Mo|sKF„-Mo| 

if +  1 

iL'  = 


it 

iL  +  \ 


(A-9) 


The  initial  phase  6^  is  estimated  from  the  requirement  that  the  total  phase  angle  at  time 
FJRST^  must  be  cither  0  or  ;r,  depending  upon  the  sign  of  the  derivative  of  the  signal. 


J  FIRST 
-(t>^FIRST 

A.1  Estimation  It«ratkHi 


forV^fr>V.j:^^  (A- 10) 


Once  the  parameters  of  the  signal  are  estimated,  we  calculate  an  estimator  of  the  sine  wave 
and  the  measurement  errors  of  each  data  point.  We  then  adjust  the  estimate  to  minimize  the 
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errors.  After  the  yth  iteration,  the  estimate  of  the  total  signal  is  Yp,  and  the  error  in  the 
estimate  is  Ep,  as  given  by 

Yp  -  Mj  +  Aj  sin|fi)yi  +  )  (A- 1 1 ) 

Ep^Yp-V^  i^lX,N.  (A-12) 

Substituting  equations  (A-1)  and  (A-2)  into  (A-12)  yields 

Ep  =Mj-M  +  AjSin^cOji  +  coi+6)-  (A- 1 3) 

Defining  AM j,  AAj,  AcOj,  and  AOj  as  the  errors  in  the  parameter  estimates  after  the  y  th 
iteration  gives 

AMj  =  Mj  -  M,  AAj  =  Aj  -  A, 

AWj  =  <0j  -  <0,  ABj  ^Bj-e.  (A- 1 4) 

Expanding  equation  (A- 13)  in  terms  of  these  errors  and  using  the  small-angle  approximation 
for  the  sine  of  the  phase  errors  results  in  the  following: 

Eji  =  AMy-i-AAysin|tt)yi  +  0yj+ Ay|Aciiyj  +  A0yjcos|a)y  +  ©yj-ej-.  (A-15) 

Equation  (A-15)  shows  tliat  the  error  in  estimating  the  fundamental  can  be  approximated  by 
an  error  in  the  offset,  plus  a  sine  wave  that  is  proportional  to  the  error  in  the  amplitude,  plus  a 
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cosine  wave  that  is  proportional  to  the  sine  of  the  error  in  the  total  phase  estimaie. 

Calculating  the  direct  current  (dc)  and  fundamental  terms  of  the  discrete  Fourier  transform  of 

Yji  allows  estimation  of  the  error  in  the  offset,  amplitude,  and  average  phase  angle.  To  avoid 

introducing  OTors  caused  by  the  partial  cycles  at  the  end  of  the  data  record,  the  Fourier 
transform  and  correlation  coefficient  arc  evaluated  only  for  the  data  points  corresponding  to 

NC  full  cycles,  i.e.,  for  iF  ^i^UL.  We  obtain  an  estimate  for  AMj,  AAj,  AQ)j,  and  Adj 
from 


(A- 16) 


(A-17) 


24  ^  /  \ 

Aaj  = - Y  I  EJi+i)cos{<t)A+eA 

^  A,NP^i=iF  ^  ^  ^ 


(A-18) 


24  iE 


24  \  - 


(A-19) 


where  the  mean  value  of  i  is  i ,  and  NP  is  the  number  of  points. 


l  =  {iF  +  iL)l2 


^/»  =  iL-iF+l 


(A-20) 

(A-21) 
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At  each  iteration,  the  estimate  of  the  signal,  and  the  error,  for  each  data  point  arc 

calculated  with  equations  (A-11)  and  (A-12).  Equations  (A'I6)  to  {A-19)  arc  used  to 
evaluate  the  entws  in  the  parameter  estimators,  which  are  then  incremented. 

This  algoridun  works  remarkably  well  for  most  digitized  waveforms.  It  rapidly  converges  to 
a  good  approximation  of  the  fundamental,  and  the  errors  approach  zero  as  j  increases.  The 
number  of  iterations  required  depends  upon  the  desired  accuracy,  the  signal  frequency,  and 
the  noise  and  distortion;  but  mainly,  it  depends  upon  the  number  of  data  points.  Computer 
simulations  indicate  that  the  results  have  converged  in  six  iterations  for  a  data  record  of  64 
points,  but  only  two  iterations  are  needed  for  a  data  record  4096  points  long. 

After  the  sine-wave  curve  fit  has  converged,  the  RMS  value  of  Ej^  determines  the  signal-to- 

noise  and  distortion  ratio  of  the  A/D  converter.  The  harmonic  distortions  can  be  computed 
via  the  Fourier  transform.  If  these  distortion  terms  are  subtracted  out,  the  remainder  is  the 
noise  of  the  converter.  If  the  noise  squared  is  correlated  with  the  square  of  the  rate  of  change 

of  signal  voltage,  i.e.,  cos^{ciM  +  6),  the  correlation  coefficient  is  proportional  to  the  total 

2 

jitter  in  the  A/D  converter  plus  the  test  equipment.  If  the  noise  is  correlated  with  (i  -  T)  ,  it 

means  that  there  is  a  small  drift  in  the  relative  frequencies  of  the  signal  and  clock  oscillator 
during  the  test. 
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APPENDIX  B 

FORTRAN  LISTING  OF  THE  PROGRAM  MONO 
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PROGRAM  MONO 
C 

C  TfflS  PROGRAM  FITS  A  CLIPPED  COSINE  WAVE  TO  64K 

MEASUREMENTS 

C 

INTEGER*2  I1(512),Y(65536).YY(65536) 

INTEGER  POINT.U  APPEAR 1(4096),APPEAR2(4096) 

REAL  V(65536)ANG2J»HI 

REAL  W{65536).SUM1(4096).SUM2(4096) 

REAL*8  ANGJDPIX)ELTA,W,THETA 
REAL*F  ONEEIGHT 

DIMENSION  A(3,3)AINV(3,3).WKAREA(18) 
CHARA'CTER*4  LAST(9) 

CHARACTER*  12  FNAME(9) 

CHARACTER*5  FILE,COMPANY 
C  INITIALIZATION  OF  STATIC  VARIABLES 
ONE=LO 
EIGHT=8.0 

DPI=EIGHT*DATAN(ONE) 

NPNTS=65536 

HNPNTS=NPNTS/2 

NC=20481 

W=DPI*NC/NPNTS 

LAST(1)='AS' 

LAST(2)=’X)AT 

LAST(3)='.64’ 

LAST(4)='.MEA' 

LAST(5)=’.CAL' 

LAST(6)='X)r 

LAST(7)='X)2' 

LAST(8)='.ONF 

LAST(9)='.ZER’ 

COMPANY^DDC 
DO  1=1,4096 
APPEAR1(I)=0 
APPEAR2(I)=0 
SUMia)=0 
SUM2a)=0 
ENDDO 
C 

OPEN(UNIT=l,nLE=T)DC_FILES.DAT,STATUS=’OLD') 

C 

C  READ  IN  DATA 
C 

5  READ(1,*END=99)FILE 
WRITE(6,*)FILE 
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INITIALIZATION  OF  DYNAMIC  VARIABLES 


C 
C 

c 

TO=0 

TS=0 

TSS=0 

TC=0 

TCC=0 

TSC=0 

RO=0 

RS=0 

RC=0 

DC=0 

S=0 

C=0 

AMP=0 

PH1=0 

U=0 

DELTA0=0 

ANG=0 

B1=0 

B2=0 

M=0 

DOI=LNPNTS 

va)=o 

wa)=o 

Ya)=o 

YYa)=0 

ENDDO 

DOJ=l,512 

n(J)=0 

ENDDO 

C 

FNAME(i)=FILE//LAST(l) 

FNAME(2)=FILE//LAST(2) 

FNAME(3)=FILE//LAST(3) 

FNAME(4)=HLE//LAST(4) 

FNAME(5)=FILE//LAST(5) 

FNAME(6)=FILE//LAST(6) 

FNAME(7)=FILE//LAST(7) 

FNAME(8)=FILE/A-AST(8) 

FNAME(9)=FILE//LAST(9) 

OPEN(lJNIT=71JFE.E=FNAME(3),STATUS=’OLD') 

C  0PEN(UNIT=73JTLE=FNAME(1),STATUS='NEW’, CARRIAGE 

CONTROL^TJST) 

0PEN(UNIT=77JTLE=FNAME(2),STATUS=’NEW’) 

OPEN(UNIT=llJTLE=FNAME(4),STATUS=’NEW) 

C  OPEN(UNIT=12,FILE=FNAME(5),STATUS=’NEW’) 
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c 

C  READ  BINARY  DATA  IN  FROM  FILE 
C  THIS  READ  FORMAT  IS  UTILIZED  SINCE  ORIGINAL  DATA 
IS  IN  DOS  FORMAT 
C 

DOK=l,128 
READ(71,’(5l2A2y)Il 
DO  1=1,512 
M=(K-1)*512+I 

Y(MH1(D 

IF  aiABSaiO))  .LT.2047)  THEN 
C  WRITE(73,*)I1(I) 

PHI=M*W 

B1=SIN(PHI) 

B2=COS(PHI) 

TO=TO+l 

TS=TS+B1 

TSS=TSS  +  B1*B1 

TC=TC+B2 

TCC=TCC  +  B2*B2 

TSC=TSC+B1*B2 

RO=RO+Il(I) 

RS=RS+Iia)*Bl 

RC=RC+I1(I)*B2 

ELSE 

ENDIF 

ENDDO 

ENDDO 

C 

C  THE  FOLLOWING  ROUTINE  INVERTS  A  MATRIX 
C 

A(l,l)=TO 

A(L2)=TS 

A(L3)=TC 

A(2,1)=TS 

A(2,2)=TSS 

A(2,3)=TSC 

A(3,1)=TC 

A(3,2)=TSC 

CALL  LINRG(3w\,3^INV,3) 

C  CALCULATE  DC,  S  AND  C 
C 

DC=AINV(  1 , 1  )*RO+AINV(  1 ,2)*RS+AINV(1 .3)*RC 
S=AINV(2,l)*RO+AINV(2,2)*RS+AINV(2,3)*RC 
C=  AINV(3.1)*RO+AINV(3,2)*RS+AINV(3,3)*RC 
C 
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C  DETERMINE  AMPLITUDE  AND  PHASE 
C 

AMP=(S**2+C**2)**0.5 
IF  (C  XT.  0)  THEN 
PHI=ATAN(S/-C)  +  DPI/2.00 
FT  SF 

PHI=ATAN(-S/C) 

ENDIF 

C 

C  FIND  COSINE  WAVE 

C 

DELTA=DPI/NPNTS 
U=PHI/DELTA 
DELTAO=PHI-U*DELTA 
IF(U.LT.  1)U=NPNTS+U 
C  WRITEdl  *)DELTAO 
DOI=l,NPNTS 
U=U+NC 

IF  (U  .GT.  NPNTS)  U=U-NPNTS 

ANG=DELTA*U+DELTAO 

Va)=DC+AMP*COS(ANG) 

YY(U)=Y(I) 

VV(U)=V(I) 

C  WRITE(12,*)Ya).V(I),U 
ENDDO 
C 

DOU=l.NPNTS 
WRITE(1 1  ♦)YY(U),  W(U) 

ENDDO 

C 

C  DETERMINE  AVERAGE,SUM  AND  OCCURRENCE  OF  EACH 
POINT 

C  REAL  RANGE  IS  FROM -2048  TO  2047.  INDEX  OF  ARRAY  IS 

TO  SERVE  AS 

C  THE  SPECDFIC  POINT  THUS  INDEX  =  POINT  +2049  => 

INDEX(1,4096) 

C  ALSO  TWO  LOOPS  ARE  SETUP  TO  DISTINGUISH  BETWEEN 

INCREASING 

C  AND  DECREASING  DATA. 

C 

D0U=1JINPNTS 
INDEX=YY(U)+2049 
APPEAR  1  (INDEX)=APPEAR1  (INDEX)+1 
SUM  1  (INDEX)=SUM1  (INDEX)+VV(U) 

ENDDO 

ISTART=HNNPNTS+1 
DO  U=ISTART,NPNTS 
INDEX=YY(U)+2049 
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APPEAR2(INDEX)=APPEAR2(INDEX)+ 1 
SUM2aNDEX)=SUM2(INDEX)+VV(U) 

ENDDO 

DATA  TO  BE  PRINTED  TO  FILE.DAT 

WRITE(77,*)'COMPANY=’.COMPANY 

WRrTE(77*yFILE=',FNAME(3) 

WRITE(77.*)TO=',TO 

WRITE(77,*)TS=',TS 

WRrnE(77*)TSS=’,TSS 

WRITE(77.*)TC='.TC 

WRITE(77,*)TCC=',TCC 

WRrnE(77.*)TSC=’.TSC 

WRrrE(77.*)'RO=’,RO 

WRrnEa7,*)'RS=’J[lS 

WRnEa7.*)'RC=’,RC 

WRrrE(77,*)'OUTPUT  OF  INVERTED  MATRDC’ 

WRITE(77,*)AINV(1 , 1 ),  AINV(1 ,2)  AINV(  1 .3) 

WRITE(77,*)AINV(2, 1  ),AINV(2,2)  AINV(2,3) 

WRITEa7,*)AINV(3, 1  ),AINV(3,2)^INV(3,3) 

WRrrE(77*)'DC=',DC 

WRITE(77,*)’S='.S 

WRITE(77,*)'C=',C 

WRnB(77*)'AMPUTUDE=’AMP 

WRITE(77,*)'PHI=',PHI 

GOTO  5 

PRINT  OUT  SUM.  APPEARANCE  AND  AVERAGE  FOR  EACH 


OPEN(UNIT=13,FILE=FNAME(6).STATUS='NEW’) 
OPEN(UNIT=14,FILE=FNAME(7),STATUS='NEW’) 
OPEN(UNIT=15,FILE=FNAME(8),STATUS='NEW’) 
OPEN(UNIT=16,FILE=FNAME(9),STATUS='NEW’) 
OPEN(UNIT=  17,FILE=’CODES  1  .DP6’,STATUS=’NEW’) 
OPEN(UNTr=18,FILE=’CODES2.DP6’,STATUS=’NEW’) 
OPEN(UNIT=19,FILE='NARROW2.DP6’,STATUS=’NEW') 
OPEN(UNIT=20JTLE='VTDE2.DP6’,STATUS=’NEW’) 
ERROR1=0 
ERROR2=0 
DO  1=1,4096 
POINT=I-2049 

IF  (APPEARia)  .NE.  0)  SUM1(I)=SUM1(I)/APPEAR1(I) 
C  WRITE(13,*)POINT,APPEARia),SUMia) 

IF  (APPEAR2a)  .NE.  0)  SUM2(I)=SUM2(I)/APPEAR2(I) 

ERRORl=SUMl(I)-POINT+ERRORl 

DIFF1=SUM1(I)-SUM1(I-1) 
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ERROR2=SUM2a)-POINT+ERROR2 

DIFF2=SUM2(I)-SUM2(I-1) 

WRITE(17,*)POINTJDIFFl 

WRnE(18,*)POINTJ)IFF2 

C  WRITEd  7,*)POINT,SUM  1  (I),  APPEAR  1  a).SUM  1  (I)-POINT 

C  WRITE(18.*)POINT,SUM2(I)^PEAR2a).SUM2(I)-POINT 
C  WRITE(14*)POINT^PEAR2a).SUM2a) 

C  IF  (APPEARia)  .EQ.  1)  WRrrE(15,*)l,POINT 

C  IF  (APPEAR2a)  .EQ.  1)  WRITE(15.*)2,POINT 

C  IF  (APPEARia)  .EQ.O)WRITE(16,*)1,POINT 

C  IF  (APPEAR2a)  .EQ.  0)  WRrrE(16,*)2,POINT 

C  IF  a  .ME.  1)  THEN 

C  CWl=SUMia)-SUMl(I-l) 

C  CW2=SUM2(I)-SUM2(I-1) 

C  IF(CW1  iT.0.1)THEN 

C  WRITE(17,*)  POINT,  APPEARia-l),SUMl(I-l), 

C  &  APPEAR1(I),SUM1(I), 

C  &  APPEAR1(I+1),SUM1(I+1)/APPEAR1(I+1) 

C  ELSE 

C  IF(CW1  .GT.  1.9)THEN 

C  WRITE(18,*)  POINT,  APPEARia*l),SUMl(M), 

C  &  APPEAR1(I),SUM1(I), 

C  &  APPEARia+l),SUMl(I+l)/APPEARl(I+l) 

C  ENDIF 

C  ENDIF 

C  IF(CW2I.T.0.1)THEN 

C  WRnE(19,*)  POINT,  APPEAR2a-l),SUM2(I-l), 

C  &  APPEAR2a),SUM2(I), 

C  &  APPEAR2(I+l),SUM2a+l)/APPEAR2(I+l) 

C  ELSE 

C  IF  (CW2  .GT.  1.9)  THEN 

C  WRITE(20,*)  POINT,  APPEAR2a-l),SUM2(I-l), 

C  &  APPEAR2(I),SUM2(I), 

C  &  APPEAR2(I+1).SUM2(I+1)/APPEAR2(I+1) 

C  ENDIF 

C  ENDIF 

ENDDO 

ERROR1=ERROR1/4096 
WRITE(17,*)'AVERAGE  ERROR=’3RRORl 
ERROR2=ERROR2/4096 
WRrrE(18,*)’AVERAGE  ERROR=’3RROR2 
CALL  EXIT 
END 


111 


APPENDIX  C 
STATISTICAL  ANALYSIS 
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C.l  INTRODUCTION 


For  the  ideal  case  of  a  noise-free  analog-to-digitaJ  (A/D)  converter,  from  N  =  512K  collected 
data  points  we  can  determine  the  code  edges  of  an  Af  bit  A/D  converter  to  a  resolution  of 
fN  =  0.03  least  significant  bit  (LSB).  To  analyze  the  data,  we  sort  the  observations  in 
order  of  increasing  input  signal  voltages,  and  then  determine  the  largest  and  smallest  input 
voltages  conesponding  to  each  possible  output  code.  With  this  data  reordering,  a  tiK)notonic 
and  noise-free  A/D  converter  would  produce  a  monotonic  output  code  sequence. 

For  an  A/D  converter  with  internal  noise,  the  problem  of  determining  the  differential 
nonlinearity  and  monotonicity  is  complicated,  as  illustrated  in  figures  C-1,  C-2,  and  C-3.  For 
purposes  of  illustration,  we  consider  an  A/D  convener  with  internal  noise,  and  we  assume 
that  one  of  the  code  edges  has  an  error.  The  nominal  voltage  of  the  lower  edge  of  code  2  is 
+1.5  mVolts.  However,  we  assume  that  code  2  has  an  error  of  -0.75  mVolt.  The  heavy  line 
in  figure  C-1  shows  the  output  code  as  a  function  of  input  signal  voltage  in  the  absence  of 
noise.  Code  1  is  assumed  to  be  narrow,  with  a  width  of  only  0.25  LSB,  while  code  2  is 
assumed  to  be  1.75  LSB  wide.  If  the  internal  noise  of  the  A/D  convener  is  assumed  to  be 
Gaussian-distributed,  the  probability  that  the  output  code,  x,  equals  or  exceeds  i  is  given  by 

(C-1) 


where 


X  =  the  output  code, 

Vin  =  the  input  voltage, 

Vi  =  the  lower  code  edge  of  code  i 
o  =  the  RMS  noise  of  the  A/D  convener,  and 
F  =  the  cumulative  normal  error  function. 


1 

VTi 


u 


dz 


(C-2) 
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Figure  C-1.  Output  Code  vs.  Input  Voltage  for  an  A/D  Converter  with  0.5  LSB  Noise 


Figure  C-2.  Probability  of  Code  vs.  Input  Voltage  for  an  A/D  Converter  with  0.5  LSB  Noise 
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Figure  C-3.  Ratio  of  Variance  to  the  Mean  of  the  Ocxurrenccs 
vs.  Ratio  of  Code  Width  to  Noise 
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The  probability  that  the  output  code  is  i  is 


P[jc  =  I]  =  PU  Si  I]  -  P[x  ^  /  +  1  ].  {C-3) 

Figure  C-2  shows  the  probability  of  outputting  various  codes  as  functions  of  input  voltage  for 
the  case  where  the  root  mean  square  (RMS)  internal  noise  of  the  converter  is  0.5  LSB.  The 
probability  density  function  of  each  code  is  a  symmetrical  bell-shaped  curve,  centered  at  the 
middle  of  the  code.  Codes  of  nominal  width,  such  as  codes  0  and  3  in  the  example,  have  a 
peak  probability  of  68  percent.  Narrow  codes  have  a  lower  probability,  while  wide  codes 

hav  a  higher  probability  and  a  flattened  peak.  From  the  probability  functions  shown  in 
figure  C-2,  the  expected  value  of  the  output  code  £[jc]  over  all  input  voltages,  Vin,  is 

calculated  from 


2047 

£[r:]=  ^iP[x  =  i].  (C-4) 

i=-2028 


For  each  input  voltage  Vin,  expected  or  (mean)  output  code  is  given  by 

2047 

£[x]=  X  /Plx  =  i  I  V  =  Vw]  (C-5) 

i=-2048 

where  Pr[x  =  i  I  V  =  Vin  ]  is  the  conditional  probability  of  occurrence  of  code  i,  given  the 
input  voltage  is  Vin.  fhe  mean  output  is  shown  in  figure  C-1  as  a  thin  line  which  is  nearly 
straight,  except  for  a  slight  bump  caused  by  the  assumed  code-edge  error. 

The  detailed  derivation  of  the  expected  value  and  variance  of  the  number  of  occurrences  of 
code  i  in  a  sample  of  size  N  of  input  voltages  Vin  is  described  in  the  following  sections 

(written  by  J.  E.  Kriegel).  The  analysis  reveals  that  the  expected  value  of  the  number  of 
occurrences  of  code  i ,  E\Hi\  is  equal  to  the  density  of  input  samples  {D  in  samples/LSB) 

times  the  code  width  (V,-4.  i  -  Vi). 
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£[W/]=D  (Vi^l-Vi), 


(C-6) 


The  variance  of  the  number  of  occurrences  will  depend  upon  the  expected  value  (of  the 
number  of  occurrences)  and  the  ratio  of  code  width  to  noise.  Let  u  denote  the  ratio  of  code 
width  to  noise  and  S(u)  denote  the  ratio  of  the  variance  to  the  mean  of  occurrences. 

Rgure  C-3  shows  a  plot  of  5(m)  versus  u. 

Rx>m  the  spectral  analysis  test  we  determined  that  the  Data  Device  (Corporation  (DDC)  A/D 
converter  has  an  RMS  random  noise  of  approximately  0.5  LSB.  Hence,  for  any  given  input 
voltage  the  output  code  can  vary  randomly  over  a  small  range.  As  a  result  of  the  random 
noise,  the  data  analysis  employed,  which  is  based  upon  statistical  methods,  will  be  less 
precise  than  for  the  noise-free  case,  and  will  have  an  associated  confidence  level  of  less  than 
100  percent.  Die  values  of  code  centers  and  edges  calculated  from  this  analysis  will  have  a 
standard  deviation  on  the  order  of  0.1  LSB.  In  our  analysis  we  are  calculating  8  K  code 
widths  (4  K  each  for  positive-  and  negative-going  signals),  half  of  which  are  only  based  on 
38  to  50  observations.  From  this  analysis  it  can  be  concluded  that  one  or  more  calculated 
values  of  code  width  may  be  in  error  by  up  to  0.4  or  0.5  LSB. 


C2  STATISTICAL  ANALYSIS:  THE  LOCKED  HISTOGRAM  TEST  METHOD 

In  the  locked  histogram  test  method  for  an  A/D  converter,  the  input  sample  voltages  are 
derived  from  a  cosine  input  and  form  a  predetermined  set  of  values  in  which  there  is  a 
constant  Cocked)  phase  difference  between  consecutive  values.  For  each  input  voltage,  the 
converter  produces  an  output  code;  and  for  an  input  set  of  voltages,  a  histogram  of  the 
number  of  occurrences  of  each  code  can  be  created.  This  histogram  depends  not  only  on  the 
range  and  number  of  input  voltages,  but  also  on  the  design  and  internal  noise  of  the 
converter. 
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C2.1  INPUT  SAMPLE 


To  evaluate  the  device,  an  input  sample  of  N  (=  2*6  say)  voltages,  is  produced  by  the 
phase-locked  system 

V,„  =  V^cos(a  + /1 2  tt/AO,  oc  constant,  n  =  0, 1, ....  N-1.  (C-7) 

A  value  of  =  2.171  volts  is  used  in  the  test.  It  is  assumed  that  input  samples  of  size  N  are 
created;  the  phase  difference  between  consecutive  values  of  V,;,  in  any  of  the  samples  is 
locked  at  2  Jt  /  N.  Without  the  phase-lock  feature,  any  sample  of  size  N  should  be  regarded 
as  arising  from  a  random  process.  One  such  process  is 

V^  =  VoCos((j»  (C-8) 

where  4>  is  uniformly  distributed  on  [0, 2  jtl- 

C  J.2  OUTPUT  AND  THE  A/D  CONVERTER 

For  each  input  voltage,  V|n.  the  A/D  convener  outputs  an  integer  value  (a  code)  in  the  range 
[-2048, 2047].  For  a  particular  fixed  voltage  repeatedly  input  to  the  convener,  the  output 
code  value  may  vary  due  to  internal  noise  which  is  characterized  by  a,  the  RMS  noise 
voltage.  The  probability  of  obtaining  an  output  code  value  x  =  i  for  a  particular  input  is 

written 

P[j:  =  t|V/„]  (C-9) 

in  which  is  to  be  regarded  as  an  index  associated  with  a  particular  value  of  an  input 
voltage.  The  accepted  model  for  equation  (C-2)  is 
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(C-10) 


P[x^i  I  Vin]  =  P[x^i  I  Vin]  -P[acS/  +l|  V,>,  ] 


1 

VTi 


Vin  -  Vi+ 1  ^in  -  V'l 

where  A  = - - - ,  B  =  — - —  and  V,-  =  the  lower  code-edge  voltage  for  code  i . 

Note  that  V,-  has  the  property  that 


P[x^i  !  Vin Vi]  =0.5. 

It  is  convenient  to  assume  in  this  analysis  that  the  initial  voltage  phase  a  is  zero.  The  value 
of  a  does  not  play  a  significant  role,  and  the  algebra  is  significantly  reduced  for  a  =  0.  With 
a  =  0,  and  introducing  co  = 


A  =  A(n  ©)  =  ^[Vq  cos(n  co)  -  Vi+ 1 } 
B  =  Bin  0))  =  ^  [Vo  cosin  to)  -  V/) 
where  V  s  Vin  =  Vocos(/ico),  n  =  0, 1, .... /V- 1. 

C3  THE  OUTPUT  CODE  HISTOGRAM 


(C-11) 


The  (expected)  value  E(i)  of  the  histogram  height  at  output  code  i  for  repetitions  of  the  input 
sample  of  voltages  is  given  by 
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£(/)=XP[^  =  fc] 


N~l 


n=0 


V2¥ 


’B 


dt. 


(C-12) 


Setting  //(nto)  =  -pLr 
\2jt 


B 


e-i^/2  which  is  a  well-behaved  function  of  n  to. 


N-\ 

E(i)=  Ydf^inta) 
n=0 

N~l 

=  i  X  H{na) 


Cl) 


n=0 

N-1 


(C-13) 


=  57  ]E//(«ci))to. 

n=0 


Rom  the  definition  of  a  definite  integral, 


f27c 

^  //(0)d0  =  /™^Z//(0r)(er+l-er) 


(C-14) 


N-  1 

in  which  0  =  0o  <  0i ...  <0^  =  271  and  max  (0^^.]  -  0,)  0  as  ->  oo,  it  follows  that,  to 

r  =  0 

significant  accuracy. 


E(i)  = 


jV 

2n 


flK 

//(0)  dQ 


(C-15) 
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since  "0,"  =  r  O),  "6^+1  -  0/'  =  O),  and  N  is  large.  Hence, 


E(i)  = 


N 

JC 


Pl  ^  r. 


0 


V2  n 


•5(6) 


.-r2/2 


>^(0) 


dr  d0 


"Vo 

r  ' 

_  a 

* 

1 

~  Tl 

•1 

-Vo 

V2  71 

V-V/+1 

^JVo^ 


MV. 


(C-16) 


Since  the  interior  integral  is  the  difference  of  two  error  functions,  the  interior  integral  is  zero 
for  V  >  Vj+j  +  ka;  and  V  <  V,  -  kc,  for  Jt  ^  5. 


Hence, 


£(0  = 


N 

*Vo 

< 

f  c 

1 

n 

« 

-Vo 

V2* 

/•  Yz^ 

o 

V-V/+1 


(Jj 


yjVo^V^ 


MV 


(C-17) 


Vi+l-^c 

r 

f  '  1 

_  ^ 

< 

1  1 

« 

Vi  -ka 

m 

V-Vi+l 

cr 

VVo2-V2 

dV; 


and  by  the  use  of  a  standard  mean-value  theorem  for  integrals  of  products  of  functions, 

V-Vj 

I 

Ki  dV 


E(l)  =  - 


Vi-k<s 


1 


V2¥ 


V-  Vi+I 


e-t^f2  dt 


(C-18) 


where 
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min  Ki  S  max  foT  Vi  -ka  V  ^  Vi+l  +kc!.  (C-19) 

V  ^Vo^V2  *  V 

Note  that  since  ~  2  volts  and  a  and  V/+]  -  V,  are  of  the  order  of  a  few  millivolts,  for  a 
given  i,  Ki  is  essentially  constant.  Equation  (C*6)  can  then  be  written  as 


E(i)  =  N 


n 


VHt 


.  v-v/ 

c 

V-V/+1 


dV. 


(C-20) 


Vi+i-Vi 

Introducing  A  = - ,  and  the  linear  substitutions 


V  ..  Vi 

2  =  1/  =  2  - 


E(i)  =  No 


r  C 

h 

1 

71 

V2¥ 

-00 

U 


.-t2/2 


dt  >  dU 


(/-A 


(C-21 


The  double  integral  in  the  square  brackets  is  a  function  of  A  only.  Omitting  the  details,  the 

K' 

integral  can  be  shown  to  have  the  value  A  by  reversing  the  order  of  integration  so  that 

E(0  =  ^  (Vi+l  -  Vi).  (C-22) 

This  gives  substance  to  the  qualitative  statement  that  the  probability  of  occurrence  of  code  i 
is  given  by  the  code  density  times  the  code  width. 
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C  J.1  THE  VARIANCE  ASSOCIATED  WITH  E(i) 


The  variance  G(0  associated  with  £(/)  is  given  by 


E(i)  =  SP[JC  =  i  I  Vin]il-P[x  =  i  I  V/„]) 
Vin 


Following  the  methods  given  in  section  C-3, 

E'i 


G(i)=  E(i)-No 


K 


/* 

oo 

- 

1 

(  V  ^  1 

d! 

2 

'  dU 

-  •> 

•oo 

V2  n 

£/-A  J 

=  E(i)-Na^I 


(C-23) 


(C-24) 


where  K'i  and  Ki  are  essentially  equal.  Note  that  the  integral  /  in  square  brackets  is  a 
function  of  A  only.  By  differentiating  I  with  respect  to  A,  it  can  be  shown  that 


£  =  where  = 


(in 


0 


(C-25) 


and  it  follows  that 


I  =  Aerf 

the  constant  term  being  obtained  by  numerical  integration  of  the  integral.  *  For  A  >10,  for 

2 

example,  it  is  found  that  I  ~  A  -  ;  hence, 

VTt 

1  R.  L.  Fante  is  gratefully  acknowleded  for  providing  an  alternative  direct  method  of 
evaluating  1. 


A^-A2/4  _  _2_ 


(C-26) 
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m 

N 


,  Vi+i-Vi 

=  a-A.  A=— - 


G(i) 

N 


/ 


=  a 


It 


A 


A>  10 


2 

=  a-7-p,  A>  10 
^  VJl 


(C-27) 


So  G(i)  is  independent  of  A  for  A  large. 


C.4  RESULTS 

The  ratio  referred  to  in  the  reference  as  S(u),  us  A,  is  given  by  S(A)  =  1  -  ^  from 
equation  (C-26). 


Behavior  of  5^  A -»  0 


Using  V  Hospital's  rule,  lim 

A  -♦0 


lim  =  0,  so  that  5(A)  ->1  as  A  0. 
A  -»0 


Behavior  of  5  as  A -»  00 
From  equation  (C-27)  for  A  >  10, 


5(A)  =  ^  — >  0  as  A  — >  . 


(C-28) 
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Behavior  of  the  "normalized  variance”  A  5(A)  as  A  -»  «> 


A  5(A)  s  A 

=  A 
2 


G(i) 

E(i) 

(1  1 


for  A  >  10 


independent  of  A  for  A  large. 


(C-29) 


The  value  of  A  5(A)  obtained  by  simulation  for  A  large  was  1.13,  which  is  the  value  of 


C J  EXTENSION  FOR  NON-PHASE-LOCKED  SAMPLES 

In  other  test  regimes,  similar  large  input  voltage  samples  are  used  in  which  the  phase  values 
(and,  hence,  phase  differences)  are  not  deterministically  specified.  Often  the  samples  are 
produced  by  the  process 


Vin  =  V^cosi^)  (C-30) 

where  <|>  is  uniformly  distributed  on  [0, 2ji].  The  above  analysis  is  easily  modified  to  include 
such  random  input  samples,  and  results  are  consistent  with  details  published  in  the  literature. 
The  additional  analysis  is  omitted  here. 


C.6  CONCLUSIONS 

This  appendix  presents  a  general  analysis  for  analyzing  the  statistics  of  narrow  and  wide 
code  widths  using  the  locked  histogram  test  method.  Estimates  of  the  mean  and  variance  of 
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the  code  widths  are  dependent  on  a  number  of  variables  which  include  input  sample  size  as 
well  as  the  RMS  value  of  the  noise  voltage  of  the  A/D  converter. 

The  above  analysis  shows  that  when  the  code  width  is  significantly  smaller  than  the  noise 
RMS  value,  the  ratio  of  the  variance  and  expected  mean  of  the  histogram  height  for  an  output 
code  approaches  1.  However,  when  the  code  width  is  much  larger  than  the  RMS  noise 

voltage,  this  ratio  behaves  like  ^code^vidS  method  of  testing,  increasing  the 

size  of  the  input  sample  does  not  change  the  results.  However,  if  the  size  sample  is  too  small 
it  will  not  produce  results  with  any  significance.  The  method  is  general  and  can  be  extended 
to  other  well-behaved  input  waveforms. 
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GLOSSARY 


ACRONYMS 

ac 

A/D 

ADC 

alternating  current 
analog-to-digital 

Analog  Devices  Corporation 

D/A 

dc 

DDC 

digital-to*analog 
din»;t  current 

Data  Device  Corporation 

FFT 

fast  Fourier  transform 

LSB 

least  significant  bit(s) 

Mwords 

megawords 

NWS 

North  Warning  System 

PC 

personal  computer 

RC 

RMS 

resistor-capacitor 
root  mean  square 

SFDR 

SINAD 

SNR 

spurious-free  dynamic  range 
signal-to-:.oise-and-distortion  ratio 
signal-to-noise  ratio 

TTL 

transistor-transistor  logic 

UAR 

Unattended  Radar 
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